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Prologue 

Random walks in random environments have been at the center of the proba- 
bilists' interest for several decades. A specific class of such random walks goes 
under the banner of the Random Conductance Model. What makes this class 
special is the fact that the corresponding Markov chains are reversible. This 
somewhat restrictive feature has the benefit of fruitful connections to other, 
seemingly unrelated fields: the random resistor networks and gradient fields. At 
the technical level, many of the problems are thus naturally embedded into the 
larger area of harmonic analysis and homogenization theory. 

This survey article is an expanded version of the set of lecture notes written 
for a course on the Random Conductance Model that the author delivered at 
the 2011 Cornell Summer School on Probability. A personal point of view pro- 
moted here is that the Random Conductance Model belongs to the collection of 
"paradigm" problems such as percolation, Ising model, exclusion process, etc, 
that are characterized by a simple definition and yet feature interesting and non- 
trivial phenomena (and, of course, pose interesting questions in mathematics). 
The text below attempts to summarize the important developments in the un- 
derstanding of the Random Conductance Model. While paying most attention 
to recent results, much of what is discussed draws on by-now classical work. 

The text retains the layout of lecture notes that have been spiced up with 
comments and references to related subjects. The general structure is as follows: 
The first section introduces the three rather different areas where the Random 
Conductance Model naturally appears. Sections 2-5 then deal predominantly 
with the first such area — namely, the various aspects of the limit behavior 
of random walks in reversible random environments. Section 6 then applies the 
introduced machinery to the remaining problems. A number of Problems are 
mentioned throughout the text; these refer to questions that are either solved 
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directly in the text or remain a subject of research interest until present day. 
Easier questions are phrased as Exercises; these are of varied difficulty but should 
all be generally accessible to graduate students. 
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1. OvervieviT and main questions 
1.1. Random conductance model 

We begin with the definition of the problem in the context of random walks in 
random environments. Consider a countable set and suppose that we are given 
a collection of numbers {uJxyjx.y^-v with the following properties: ujxy > with 

7r„(x) := uj.xy e (0, cx)), x^Y, (1.1) 

and the symmetry condition 

iOxy^ujyx, x,y&y. (1.2) 

We will predominantly take Y to be the hypercubic lattice Z'' naturally embed- 
ded in W^. The quantity oj^y is called the conductance of the pair (x, y) — the 
use of the term will be clarified in the subsection dealing with resistor networks. 

When "f has an unoriented-graph structure with edge set rf, we often enforce 
^xy — whenever (x, y) ^ in that case we speak of the nearest-neighbor 
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model. Such a model is then called uniformly elliptic if there is a G (0, 1) for 
which ^ 

a<ujxy<-, {x,y)€S. (1.3) 

When y := Z'', we use the phrase "nearest-neighbor model" for the situation 
when S' is the set of pairs of vertices that are at the Euclidean distance one 
from each other. 

The aforementioned "random walk" in environment oj is technically a discrete- 
time Markov chain with state-space Y and transition kernel 

P^{x,y):=^^, x,y&r. (1.4) 

In plain words, the "walk" at site x chooses its next position y proportionally 
to the value of the conductance oJxy The non-degeneracy condition (1.1) guar- 
antees that this chain is well defined everywhere; when positivity of tt^ fails at 
some vertices — as, e.g., for the simple random walk on the supercritical per- 
colation cluster, cf Fig. 1.1 — one simply restricts the chain to the subset of Y 
where t:uj{x) > 0. 

A key consequence of the symmetry condition (1.2) is: 

Lemma 1.1. tt^^ is a stationary and reversible measure for the Markov chain. 

Proof. Invoking the above definitions we get 

'^i^{x)Puj{x,y) = UJxy = t^yx = TT^ (j/) (y , x) , (1.5) 

which is the condition of reversibility (a.k.a. the detailed balance condition). 
The fact that tt^^ is stationary follows by summing the extreme ends of this 
equality on x. □ 

Note that for the nearest-neighbor model on 'L'^ with conductances ijjxy ~ 1 
if ja; — y| = 1 and = otherwise, the above Markov chain reduces to the or- 
dinary simple (symmetric) random walk. In this case the increments of the walk 
are i.i.d. which permits derivation of many deep conclusions — e.g., Donsker's 
Invariance Principle, Law of Iterated Logarithm, etc. However, when ui is non- 
constant, the increments of the chain are no longer independent; worse yet, they 
are not even stationary. As we will see, this can be overcome but only at the cost 
of taking w to be a sample from a shift-invariant distribution. This reasoning 
underpins the large area of random walks in random environment of which the 
above chain is only a rather specific example. 

Let be the space of all configurations {iOxy) of the conductances. This 
space is naturally endowed with a product cr-algebra ^ . A shift by x is the map 
Ta; : ri — >■ acting so that 

iTxUj)yz ■■= UJy+x.z+x, X, y, ZeZ"*. (1.6) 

We will henceforth assume that P is a probability measure on (f2, ^) which is 
translation invariant in the sense that 



PoT^i^p, xeZ'^. 



(1.7) 
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Fig 1.1. Naturally included in the family of Random Conductance Models is the simple ran- 
dom walk on a supercritical percolation cluster (for details on the definition and properties of 
percolation, see the monograph of Grimmett [72]). Here the conductances are nearest-neighbor 
only and take values either zero or one independently at random with the same probabilities 
everywhere. The density of conductance- one edges exceeds the percolation threshold so there 
is an infinite connected component of vertices joined by conductance- one edges; only the ver- 
tices in this component are retained in the figure. At each time the random walk chooses a 
neighbor of its current position at random and passes to it. The marked vertices depict those 
visited by a sample path of such random walk started at the center before it exits the box at 
the point on the right indicated by the arrow. 



We recall that this measure is said to be ergodic if F{A) G {0, 1} for any event A 
with the property t~^{A) = A for all a; € Z'^. A canonical example of an 
ergodic P would be the nearest-neighbor model where the values of conductances 
are chosen independently at random from the same distribution. We will use E 
to denote expectation with respect to P. 

Let us now turn to the main questions one may wish to ask concerning the 
above setup. For this let X = (X„) denote a sample path of the above Markov 
chain and let denote the law of X subject to the initial condition 

P^{Xo^x):=l. (1.8) 

Let denote the n-th power of the transition kernel P^^, i.e., 

PZ{^,y) = P:iXn = y). (1.9) 

The aforementioned connection with the special case of simple symmetric ran- 
dom walk leads to the following questions: 
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Problem 1.2. Does the limit 

lim — (1.10) 

n— >oo 77, 

exist almost surely? Under what conditions is it zero (as it is for the simple 
random walk)? 

Problem 1.3. Under what conditions does the path obey an invariance principle 
— i.e., does its law tend to Brownian motion under diffusive scaling of space 
and time? And if so, what is the rate of convergence? 

Problem 1.4. Does one have a local CLT as well in the sense that 

(1-11) 

whenever y can be "comfortably" reached by the random walk from x in n steps? 

As it turns out, there are subtle but important differences in the precise 
technical sense in which these asymptotic statements might be true, or at least 
provably true. Indeed, there are two natural laws on the path space that are 
considered in the literature: the aforementioned quenched law P^{—) and the 
annealed or, more accurately, averaged law EqP^{—) where Q is a specific (nat- 
ural) measure on environments (similar to P) . An advantage of the annealed law 
is that, thanks to averaging, it allows for an easier control of the irregularities of 
the environment; a drawback is that the path law under it is no longer Marko- 
vian. As we will see, one of the main challenges for the Random Conductance 
Model that prevail to the present day is the resolution of: 

Problem 1.5. Does the annealed invariance principle imply the quenched in- 
variance principle? (Here and henceforth the words annealed and quenched des- 
ignate the path distribution that is considered for the scaling limit.) 

We remark that, for general random walks in random environments, the an- 
nealed and quenched law can be dramatically different. See Fig. 1.2. 

As soon as the above "fundamental" questions have been resolved, one can 
try to imitate various derivations that have over years been accomplished in 
the context of the simple random walk. This leads to further rather interesting 
questions, for instance: 

Problem 1.6. What are the intersection exponents — i.e., the decay exponents 
for the probability of non-interection up to the first-exit time from a ball of a 
large radius — of several independent copies of such random walks ? 

Problem 1.7. Does the (chronological) loop-erasure of the walk have the same 
scaling limit as the simple random walk? And how many steps of the walk are 
needed to generate n steps of the loop- erased walk? 

Problem 1.8. Is there a scaling limit for the trace of the walk in Fig. 1.1 as 
the size of the box tends to infinity? 

The last question naturally puts us into a bounded domain where, as it turns 
out, many additional technical difficulties arise compared to the full lattice. 
However, even the following questions are quite relevant: 
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Fig 1.2. j4n example of a random walk in a random environment where, at each vertex, one 
of the North or East arrows is chosen independently at random (with equal probabilities). 
The random walk is then forced to follow the arrows. The extremity of this example is seen 
from the fact that while the quenched law of the path is deterministic — and no invariance 
principle can hold for fluctuations — the averaged law looks like an ordinary North & East 
random walk whose fluctuations are described by the Central Limit Theorem. 

Problem 1.9. Is there a scaling limit of the random walk among random con- 
ductances restricted to the half-space, quarter .space or a wedge (i.e., for the 
problem with conductances "leading" outside these regions set to zero)? 

1.2. Digression on continuous time 

Although the discrete-time Markov chain is very natural, one is often interested 
in a continuous-time version thereof. We will therefore introduce these objects 
right away and discuss some of the technical issues that come up in this context. 

There arc two natural ways how to make the time flow continuously. First, 
we may simply Poissonize the discrete time and consider the transition kernel 

Ql{x,y):=Y,^^^-'PZ{^>y)- (1-12) 

ri>0 

The corresponding (continuous-time) Markov process is then referred to as 
constant- speed random walk among random conductances (CSRW), where the 
adjective highlights the fact that the jumps happen at the same rate regardless 
of the current position. 

Another natural way how to make time flow continuously is by attaching a 
clock to each pair {x,y) that rings after exponential waiting times with expec- 
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tation l/ujxy This can just as well be done by prescribing the generator 

(U/)(x):-^c.,,[/(y)-/(a:)], (1.13) 

V 

and demanding that the corresponding transition kernel is the (unique) 
stochastic solution of the backward Kolmogorov equations, 

^RL(x,y) = ^U(a;,z)RL(z,2/) (1.14) 

Z 

with initial condition 

Here 5x {z) equals one when x = z and zero otherwise. This leads to the variable 
speed random walk among random conductances (VSRW), because the resulting 
Markov chain at x makes a new jump at rate 7r^(a;). 

A specific problem with the VSRW is that the walk may escape to infinity 
in finite time — a blow-up occurs. (This will not happen for the discrete-time 
walk and thus also the CSRW.) A simple criterion to check is: 

Exercise 1.10. Consider a configuration u of conductances such that Huix) G 
(0,oo) for each x. Let {X^) be the path of the discrete-time random walk among 
conductances ui and let To,Ti, . . . be the times between the successive jumps of 
the corresponding VSRW. Show that 

The upshot of this Exercise is that the question of blow-ups in VSRW can 
be resolved purely in the context of the discrete-time walk. We refer to, e.g., 
Liggett [95, Chapter 2] for a thorough discussion of such situations. See also 
Exercise 2.8 in Sect. 2.2. 

The above transition kernels are distinguished by their invariant measures 
and natural function spaces they act on. Indeed, we can write R^ as 

Rl{x,y) := ((5y,e*'-"4)^,(2,^, (1.17) 

where we think of £^(Z'^) as endowed by the counting measure. On the other 
hand, the constant speed Markov chain admits the representation 

Qli^,y) ■■= ^(<5„e*(P"-i)4),.(^^) (1.18) 

where 1'^{'k^) is the space of functions /: ^ M that are square integrable 
with respect to the measure tt^j on . In this case the generator of the Markov 
chain is simply P^j — 1. The reason why one uses different underlying measure 
on 1/^ in the two cases is seen via: 
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Exercise 1.11. Show that is symmetric oni'^{'L'^) while Pcj — 1 is symmetric 
on Piii^). In particular, the VSRW is reversible with respect to the counting 
measure on Z'* while the CSRW is reversible with respect to tt^ . 

It is clear that the constant-speed chain will follow the discrete-time chain 
very closely, but the variable-speed chain may deviate considerably because 
its time parametrization depends on the entire path. This discrepancy will be 
particularly obvious in the places where, in comparison with the neighbors, 
TTuiix) is either very small (VSRW gets stuck but CSRW departs easily) or very 
large (VSRW departs easily but CSRW gets stuck). This may or may not be a 
disadvantage depending on the context. 

1.3. Harmonic analysis and resistor networks 

The above (discrete-time) Markov chain is in a class of models for which we 
can apply a well-known connection between reversible Markov processes and 
harmonic analysis/electrostatic theory. This connection goes back to the work of 
Kirchhoff in mid 1800s (KirchhofF [87]) and it underlies many modern treatments 
of Markov processes. For our purposes the best general introductory text seems 
to be the monograph by Doyle and Snell [48]. 

We begin by introducing some relevant notions for the full lattice; the finite- 
volume counterparts will be dealt with later. For a configuration of the conduc- 
tances (uJxy) and a function /: Z'' R, let us define 

f(/):=l^c..,[/(y)-/(x)]l (1.19) 

In physics vernacular, this is the electrostatic or Dirichlct energy corresponding 
to the electrostatic potential /. We then define the effective (point-to-point) 
resistance R{x,y) between x and y by the formula 

R{x,y)-' := inf{f (/) : / € £\n^), f{x) = 1, f{y) = O}. (1.20) 

More generally, we define an effective point-to-set resistance R{x, A) by requiring 
f{y) = for all y € A in the formula above. Of course, both £{f ) and R{x, y) 
depend on w, but we leave that notationally implicit. 

A key problem now is a computation, an analysis of various scaling proper- 
ties, of the effective resistance. As a warm-up, consider now the homogeneous 
problem when the conductances are equal to one for nearest neighbors and zero 
otherwise. Leaving aside some technical issues, any minimizer of the Dirichlct 
energy in (1.20) will then obey 

E [/(")-/W]-0' vel^Mx^y}, (1.21) 

u: \u—v\ — l 



with 



/(y) - fix) = 1. 



(1.22) 
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Fig 1.3. An example of an electrostatic problem connected to the Random Conductance Model. 
Here part of the percolation cluster in a slab with vertical coordinates in the interval [—N, N] 
is attached to metal plates with a given voltage difference. The edges present in the cluster 
have resistivity one, the edges that are absent are total insulators. A key question is to find the 
total current density — per unit area of the plates — running through the system. Another 
question is the value of the electrostatic potential at the origin. 

In other words, / is discrete harmonic everywhere away from x and y. It is an 
interesting exercise in upper-division analysis to solve: 

Exercise 1.12. Fix I e M. For the homogeneous nearest-neighbor problem, use 
Fourier transform to solve the equation 

E [f{u)~ f{v)\^l[5^{v)~6y{v)\, v&%\ (1.23) 

u: \u—v\ — \ 

and then adjust I so that f{x) — f{y) — 1. Use this to derive an integral formula 
for R{x,y). 

We can thus check that while the following problem may appear hard, it is 
at least not ill posed: 

Exercise 1.13. For the homogeneous nearest-neighbor problem on 1? , show 
without relying on Fourier transform that R{x, y) — V2 whenever x and y are 
nearest neighbors. 

Returning to the full-fledged Random Conductance Model, let us now dis- 
cuss the (somewhat degenerate) example of the supercritical percolation cluster 
depicted in Fig. 1.3. Assuming the potential is fixed to = — 1 at a conducting 
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plate at "height" — and to Lp = +1 at the corresponding plate at "height" 
+N , the question is what is the electrostatic potential right at the center. As 
before, this potential is a minimizer of the Dirichlet energy £{f) in (1.19) sub- 
ject to the conditions that f{x) := 1 when e2 ■ x > N and f{x) := —1 when 
62 ■ X < —N. Here e^ is the coordinate unit vector in the i-th lattice direction. 

What makes this problem relevant for probabilists is the existence of a direct 
probabilistic "solution:" Let t^^'' be the first hitting time of the upper, resp., 
lower metal plate, 

:= inf{n > 0: X„ • 62 = ±iV}. (1.24) 

Then the electric potential at vertex x turns out to be given by the formula 

^ix) := P:(rf ) < r!-') - P^ir^^^ < ), (1.25) 

where is our notation for the law on paths (X„) of the random walk on 
environment u such that P^{Xo = x) = 1. The key point is that the function 
(f defined by (1.25) is harmonic with respect to the generator of the continuous 
time Markov chain (1.13) with the boundary values given as above. Here a 
function is said to be harmonic at x when Li^f{x) = 0. 

Exercise 1.14. Prove the formula (1.25) by showing that such a harmonic 
function is uniquely determined by its boundary data. 

Notice that, as soon as the conductances are non-constant, there is no reason 
why the potential ip at the symmetry point should be equal to zero — as it would 
be, thanks to symmetry considerations, for the case of homogeneous networks. 
Obviously, this is quite related to Problem 1.9. 

The concept of effective resistance is closely related to the question of recur- 
rence and transience of the corresponding Markov chain. Let 

:= inf{77, > 1 : X„ = x} (1.26) 

and, for AcZ'', 

TA:=ini{n>0: X„e A}. (1.27) 

Set Aat ;= [-N, N]'^nZ'^. The chain wiU then be recurrent if P°(fo < ta^^) 1 
as iV 00 and transient otherwise. The connection with effective resistance 
shows that the tendency to recurrence decreases with increasing conductances. 
Explicitly, we have: 

Exercise 1.15. Show that the function 

, , I P^(fo < TAc ), ifx^O, „„^ 

^(a;):=i/^ .„ ' (1-28) 

is the unique minimizer of the Dirichlet energy for the boundary conditions 
corresponding to point-to-set resistance R{x,A'j^) and use this to derive 

R{x,A%r^ = ^.(0)P°(fo > rAc„). (1.29) 

Conclude that 7rij(0)P^(fo > '''A'^) *s monotone increasing in each uj^y 
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The upshot of this observation is that if uj^y < w^.y for all pairs x, y, then 

^c.(0)P^"(fo > TA^„) < ^.'(0)P^",(fo > TAo^). (1.30) 

In particular, if the random walk is recurrent in the environment uj' then so it 
is in uj, and vice versa for the question of transience. For (say) nearest-neighbor 
Random Conductance Models subject to the cUipticity condition (1.3), recur- 
rence is thus equivalent to the recurrence of the simple random walk. However, 
as soon as ellipticity is violated, interesting problems arise. 

Consider for illustration the random walk on the supercritical percolation 
cluster. There the conductances are bounded above but not below. This still 
permits us to conclude that the random walk is is recurrent in spatial dimension 
d = 2, and if it is transient in dimension d — 3, then it is transient in all 
dimensions d > 3. A key question to resolve is thus: 

Problem 1.16. Is the random walk on almost every realization of the three- 
dimensional supercritical percolation cluster transient? 

The following question should ideally be solved before tackling Problem 1.16: 

Problem 1.17. Let ujb £ {0, 1} and let '^oo(w) denote the set of vertices in 
that lie in an infinite self- avoiding path using only edges withujh = 1. Let to' differ 
from UJ in a finite number of coordinates so that lu'j^ > luj, for all b. Assuming 
that £ ^^00(1^); show that 

P^{X is transient) = P°,{X is transient) (1.31) 

and conclude that {X is transient on "^oo} is a tail event. (In particular, for 
Bernoulli lo^ 's, it is also a zero-one event.) 

There are a good number of variations on the problem depicted in Fig. 1.3, 
but here is one that has been particularly perplexing for a number of years 
— in spite of an existing solution claimed in the book of Jikov, Kozlov and 
Oleinik [83]. The formulation goes back to Kesten's monograph on percolation 
(Kesten [85]). Consider the square box Ajv := [-N,N]'^ n and let ^at be 
the set of those edges whose both endpoints lie in the infinite bond-percolation 
cluster and also in A at. Define the effective resistance 

:= inf I [/O^) - /(y)] ' : /(^) = -ei • X when ei • x = ±iv| (1.32) 

corresponding to the boundary conditions —TV on the "left" side of the box and 
+ on the "right" side of the boundary; no boundary condition is prescribed at 
the remaining portions of the boundary. It is not hard to convince oneself that 
is at most of order A^''. but identifying a precise rate is far more challenging: 

Problem 1.18. Prove that for almost every realization of the supercritical per- 
colation cluster, the limit 

lim ^ (1.33) 

exists and is independent of the realization. Characterize its value. 
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Of course, once this has been settled, one may want to go beyond a LLN-type 
of information and study the fluctuations. Interestingly, as observed already a 
while ago by Wehr [134], the variance of RJ^ is order at most N'^ — at least in 
the elliptic setting — which suggests the following question: 

Problem 1.19. Show that the law of N^'^^^[R^^ — Ei?^^^] tends to Gaussian 
as N ^ OD. 

Recently, thanks to the work of Gloria and Otto [67] , we even know that the 
variance is actually of order A^** (at least in d > 3) so the time seems ripe for 
resolving this problem as well. 

1.4- Gradient models 

The third and somewhat unexpected context in which one naturally encounters 
the Random Conductance Model is that of gradient fields. In our formulation, 
a gradient field is a collection of M-valued random variables 4>x indexed by the 
vertices x E Z"^. Wc impose the following law: 

4{d^) := ^expj - ^i^- - -^y)} n d-^- n ^4>M<I'-)- (1-34) 

^ {x,y)eB{A) ^ xeA x^A 

Here A C Z'' is a finite set and B(A) is the set of all edges with at least one 
endpoint in A. The function F: R — >■ K is the potential which we take to be 
a continuous, even function with sufficient (e.g., quadratic) growth at infinity. 
The measure depends on the values immediately outside A which are set to the 
boundary condition (j) by the product of delta-masses. 

Gradient models are ubiquitous in physical sciences where they arise as 
effective-interface models, with (p^ giving the height of a surface above a ref- 
erence plane, or in descriptions of the fluctuation fields in critical statistical 
mechanical (spin) models. A higher-dimensional variant, particularly, (j)^ G R'', 
has the interpretation of a deformation field representing the displacements of 
atoms in a crystal from their ideal positions. Further applications can be found 
in field theory and material physics. The reviews by Giacomin [64], Velenik [132], 
Funaki [55] and Sheffield [122] give more information and further connections. 

We will actually consider the measure (1-34) to be a law on the sigma- field 
of gradient events 

^:=a((/), -0j,:x,yeZ'^), (1.35) 

which is legitimate since the corresponding restriction of fj,'^ does not depend 
on the values of (j) but only on their differences. This restriction is dictated 
by practical reasons — the actual "height" of an interface is usually of lesser 
importance than the "shape" of its configuration — but also due to technical 
restrictions in low spatial dimensions. We say that a measure /i on (R^ , JF) is 
a gradient Gibbs measure (GGM) if for every A G and any finite A C 'Z'^, 



(1.36) 
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where the expectation is over the boundary condition (j). Put another way, this 
says that the conditional probabihty of A given the configuration </> outside A 
is exactly the measure (1.34). 

Before we start discussing the relevant problems arising in this subject area, it 
is interesting to note two special instances of the above formalism. The first one is 
the d = 1 case. Let us assume that A is connected and, in fact, A := {—N, N)C\'L. 
Then the law of of the gradients, {4>x+i — 4'x)x=-n ^® i-i-d. — with marginal law 
proportional to e^^^^'-'dry — conditional on 

N-l 

{(fix+l - (l>x) = <I>N - <I>-N- (1-37) 

x=-N 

This situation can be analyzed with the help of standard methods of large- 
deviation theory (cf, e.g., Dembo and Zeitouni [45], den Hollander [79]) — in 
fact, Cramer's theorem more or less suffices — and so one can prove: 

Exercise 1.20. Suppose d ~ 1 and a linear boundary condition, i.e., (px '■= 

tx for some i G M. Show that, for any continuous, even potential V growing 
superlinearly at infinity, the law of 

t^N-'/^[cj)^tNi-tN], -l<t<l, (1.38) 

linearly interpolated into a continuous function, scales to a Brownian bridge 
as N ^ oo. Characterize the variance at t = 0. 

Another instance of special interest is that when V is quadratic, 

V{v) ■■= (1.39) 

for some stiffness k > 0. In this case the above measure is Gaussian and so it 
is amenable to explicit calculations. In fact, for (say) zero boundary condition 
= 0, one can even pass to the limit A f Z'', provided one restricts to the 
sigma-algebra of gradient events (1.35). This restriction is necessary because in 
dimensions d= 1,2, the law of (pQ is not tight in this limit. To see this in more 
explicit terms, note that 

Cov^o (0^,0j^) = GA(x,y), (1.40) 

where Q/^{x,y) is the Green's function associated with the discrete Laplacian 
with Dirichlet boundary condition on dK. In probabilist's terms, GA(a;, y) is the 
expected number of visits to y by the simple random walk started at x before 
it exits from A. The classical formula 

GA(^,y) = /"^y;^:^ y (1.41) 

see, e.g., Spitzer [126] or Lawler [93], using the notation (1.26-1.27), provides 
an explicit connection to the issues discussed in the previous subsection. 

An analogue of Exercise 1.20 in d > 2 will then be: 
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Exercise 1.21. Consider the Gaussian gradient model with (1.39) with k :— 1. 
For a sample of the field {(px) from the infinite-volume limit /i := lini^^x"* Ma' 
and a smooth f: M.'^ M. with compact support and J f{x)dx ~ define 

<^,(/):=eWy' ^^^j/(,^)d3;. (1.42) 
Show that in the limit e J, 0, the law of 4>e{f) is a Gaussian Af{0, ctj), where 

a2:=(/,-A-V)L2(R.)- J dxdyf{x)G{x-y)f{y), (1.43) 

where G :~ liixLA^xd Ga is the infinite- volume Green's function. (The expression 
on the right is well-define because f e Dom(A~^).^ 

The problem is meaningful in all c? > 1 but only in d = 1 we have a hope 
to describe the limit as a (real-valued) process. This is because the limiting 
continuum object, the Gaussian Free Field (GFF), is very rough in d > 2 and, 
in fact, can only be interpreted in the sense of distribution theory — hence our 
formulation using a linear functional cj)^ in (1.42). We refer to, e.g., Sheffield [123] 
for more information on the tightness issues and other aspects of the GFF. 

Having dealt with these instructive examples, let us move on to general po- 
tentials V . A remarkable feature of gradient models is that much of what has 
already been said about the quadratic case applies to any gradient model for 
which V is uniformly strictly convex — i.e., when V is C^ with V" positive 
and uniformly bounded away from zero and infinity. (We will expound on the 
specifics in the discussion of dynamical environments in Section 4.4). Unfortu- 
nately, convex potentials are not what one typically finds in models coming from 
realistic systems and / or applications and so the last decade has witnessed a ma- 
jor push to obtain a similar level of control also for non-convex interactions. This 
has so far succeeded only partially because most of the existing techniques fail 
as soon as V is non-convex anywhere, regardless how unlikely (or energetically 
unfavorable) a configuration for which this happens may be. 

Notwithstanding, there is a family of models with non-convex V that can 
be studied by way of a connection to the Random Conductance Model. These 
models are defined generally by requiring that V be given by 

e-^t") [ p(dK)e-^'"'', (1.44) 
"'(0,00) 

where p is a positive measure on positive reals. Notice that when p is supported 
at a single point, then V is quadratic, but as soon as p has at least two points 
in its support, V can be non-convex, see Fig. 1.4. (Nontheless r] 1— V{7]) will 
always be increasing on positive ?/'s.) 

An essential feature of the assumption (1.44) is that it permits us to consider 
^^(d(/)) in (1.34) as the (/)-marginal of the measure /x^(d^dK) on M^"* x (0, oo)"(^) 
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< K2 ^ 'SI the potential is not convex. 



which is given by 

^,i(d^d^)■.= ^e-"-^'f''^^ l[d^,l[S^^{d^,) n Pi<^^-y)^ (1-45) 

xeA x^A {x,y)eB{A) 

where 

HA{(j>,n) ■.= ^ K,y{^, - (f>yf . (1.46) 

{x,v)ea{A) 

To see why this holds, introduce a "private" variable K^y = Hyx for each edge 
(x, y) £ B(A) and use the additive structure of the interaction to write the 
exponential weight in (1.34) as the exponential weight in (1.45) integrated over 
the product of the p's. A key point is that, by regarding the k^j^'s as genuine 
random variables and conditioning on their values, the law of the 0's is again 
Gaussian, albeit now with a spatially inhomogeneous eovarianee structure. 

The above constructions can be performed in infinite volume; see Biskup and 
Spohn [19] for details. We will only communicate the salient conclusions: First, 
one can represent every gradient measure fj, for the potential V in (1.44) as the 
(/)-marginal of an extended measure v on pairs of configurations (0, k) such that 
the following holds: 

(1) Conditional on the (/)'s, the individual k's are independent with n^y having 
the marginal law proportional to e~^'^"=''^'^="~'^f p{dKxy)- 

(2) Conditional on the k's, the 0's are then Gaussian with eovarianee given 
by the inverse of (the negative of) the generator 

Uf{x):^ J2 ^=cy[fiy)-fix)] (1.47) 

y: \y-x\ = l 

of the Random Conductance Model with nearest-neighbor conductances 
(i^xy)- (The mean can be characterized too, but we will discuss this in the 
proof of Theorem 6.7.) 
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(3) The K-marginal is generally strongly correlated, but if the initial gradi- 
ent measure is ergodic with respect to translations, then the extended is 
ergodic as well. 

For those familiar with the Random Cluster Model (see, e.g., the monograph 
by Grimmett [73] ) and the Fortuin-Kasteleyn represenation of the Potts model 
(Fortuin and Kasteleyn [60]), the above should be quite reminiscent of the so 
called Edwards-Sokal coupling of these two processes (Edwards and Sokal [53]). 

The structure described above offers the possibility to study the gradient 
model with non-convex interaction of the type (1-44) by conditioning on the 
k's. The proof of scaling of the gradient field to the Gaussian Free Field at large 
scales then boils down to solving: 

Problem 1.22. Let {4>x) be a collection of Gaussian fields with mean zero and 
covariance given for any g: Z'' — > K. with finite support and g{x) ^ by 



where (—Li^)^^{x,y) — the inverse of the operator — — can equivalently be 
described as the full-lattice Green's function of the random walk among nearest- 
neighbor random conductances k. Show that, for any ergodic law P on the k's, 
the random functional 0e(/) in (1.42) tends to a Gaussian random variable 
P-a.s. Characterize its variance. 

As we will see this will become even more interesting once we start discussing 
gradient fields with non-vanishing tilt. Naturally, once these basic convergence 
issues are settled one can turn to more subtle questions such as, for instance: 

Problem 1.23. For a Gaussian field with covariance (— L„)^^ with zero values 
on the boundary of a cubic domain Ajv := [^N, N]'^r\'E'^ , what is the distribution 
o/max^jgAjv 4'x? What is the scaling limit of the level sets? And how about the 
Hausdorff dimension of various exceptional sets (e.g., the so called thick points)? 

These problems have recently been studied for the homogeneous lattice GFF, 
e.g., by Bolthausen, Deuschel and Zeitouni [20], Daviaud [37], Hu, Miller and 
Peres [81], Schramm and Sheffield [119], and also for the uniformly convex in- 
teractions (Miller [104]). 

1.5. Outlook 

The upshot of the above overview that all of these problems, although quite 
varied in nature, can be reduced to specific properties of the Random Conduc- 
tance Model. In particular, many of the solutions boil down to similar technical 
questions. In the rest of these notes we will attempt to explain the main ideas 
underlying the existing solutions and point out the obstacles that are known of 
for the problems that remain unresolved. 




(1.48) 



Random Conductance Model 



311 



2. Limit laws for the RCM 

The goal of this section is to exhibit the main techniques that will allow us 
to establish the validity of the SLLN (Problem 1.2) and the Functional CLT 
(Problem 1.3) for rather general Random Conductance Models. We will take a 
very pedagogical approach that starts off by addressing the simplest non-trivial 
cases of interest while isolating, as clearly as possible, various technical issues 
that come up along the way. 

2. 1 . Point of view of the particle 

A first basic problem that arises in analyzing the Markov chain (X„) for a 
fixed realization of the environment uj is that the increments of this chain are 
not stationary. A way to mend this is to invoke the first fundamental idea 
encountered in the theory of random walks in random environment: the point 
of view of the particle. Namely, instead of making a random walk run through a 
fixed environment, we will shift the environment around so that the walk remains 
always at the origin. Technically, this amounts to representing the sequence 
(tx„w) as a trajectory of a Markov chain on the space of all environments. 

Lemma 2.1. Suppose P is translation invariant. Then (ta'„i^) is a sample from 
a Markov chain on the space of conductances 17 with the transition kernel 

^(w,dw') := J2 Pu{0,x)Sr^Udio'). (2.1) 

X 

Moreover, whenever Z := ¥,t:^{0) < oo, this chain has the stationary and re- 
versible measure 

QidLu) M!)p(d,.). (2.2) 

Proof. The fact that the kernel 3^ generates the Markov chain (rx„i^) is a trivial 
calculation. For the second part, we need to invoke a bit of L^-calculus. For any 
two bounded measurable functions / = /(t^) and g = g{uj), define 

(/,.9> :=i?Q(/M5H)- (2.3) 

This is a natural inner product in L^(Q). To show reversibility (and thus sta- 
tionarity) of Q, it suffices to show that ( /, ^^g) ~ {^f, g) for any such bounded 
non-negative /, g — in fact, indicators of measurable events would be enough. 
For that case we compute 

( /, ^5) = I ^ e(^„(0)/(w)P„(0, x)g o r,(L.)) 

1 ' / (2.4) 
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where all sums arc meaningful by positivity of all terms and the assumption that 
X^a: '^o,x is integrable. Now apply t^x under the expectation to write this as 

{f,^9) = ^Y. (r-.w)o,. . (2.5) 

X 

A key property of the environment is its symmetry (1.2) whereby we get 

{t^x^)o.x = ^-x,o = Wo -s. (2.6) 
Relabeling —x for x, we thus conclude 

{L^g) ^^Y.^[f{TxUj)Loo,xg{Lo)) (2.7) 

X 

which is, rolling back the first rewrite, exactly {S^f^g). □ 
It is not hard to check that, for any bounded /, 

( /, (id - ^).9) = ^e(c.o,. [IM - I{u:)) (.g(r,c.) - g{u))). (2.8) 

X 

This will help us solve: 

Exercise 2.2. Show that, whenever Z := EtTuj{0) < oo, the operator I£ := 
^ — id with domain 

Dom(if ) := |/ e : ^(^"'- {!{-^xOj) - f{uj)f) < ooj (2.9) 

is self-adjoint and negative semi- definite. 

Notice that the stationary measure Q and the a priori law P are mutually 
absolutely continuous; we in fact even have a very explicit expression for Q. In 
the studies of general (non-reversible) random walks in random environments 
it is (usually) not too hard to infer the existence of a stationary measure but 
a key obstacle is the absolute continuity of P with respect to Q — which we 
often need to conclude that events that occur Q-a.s. also occur P-a.s. But even 
in such cases it is unusual to have any sort of explicit handle of Q. 

These considerations move us to the question under what conditions is the 
Markov chain (tx^ijj) ergodic. In order to explain this a bit better, recall that 
a stationary Markov chain (Z^) on a general state space with stationary mea- 
sure TT can always be embedded into a Markov shift as follows: Sample the initial 
state Zq from tt and then use the Markov kernel to sample a whole forward tra- 
jectory (Zn)n>i- If need be, also use the reversed chain to sample the entire 
backward trajectory {Zn)n<o- This defines — through the Kolmogorov Exten- 
sion Theorem — a law /i on trajectories of the Markov chain. The canonical 
shift — simply use Z„ for the value of Zn^i for all n — then defines a measure 
preserving transformation. 
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This construction and the Birkhoff-Khinchinc Ergodic Theorem imply that, 
for TT-almost every Zq and almost every path of the Markov chain — in short, 
for /i-almost every trajectory — the limit 

lim - V f{Zk) (2.10) 

n— foo Tl ^ — ^ 

fc=0 

exists and is finite for any function / such that / e L^(7r). However, we often 
wonder whether this limit is in fact almost surely constant — and this will only 
be true for a general / if the chain is ergodic. Explicitly, the above Markov chain 
is ergodic if any measurable set of trajectories A satisfies /i(A) G {0, 1}. 

Ergodicity will in our context be guaranteed by the following condition: 

Proposition 2.3. Suppose P satisfies the following conditions: 

(1) P(7ri^(0) > 0) = 1 an(iE7r„(0) < oo — i.e., Q exists and is equivalent to P. 

(2) V is irreducible in the sense that, for every x € Z'*, 

P( w: supP^(0,a;) > ) = 1. (2.11) 

V n>0 / 

(3) P is ergodic with respect to the translations of ifi — i.e., P(yl) G {0,1} 
for any event A such that Tx{A) ~ A for all x. 

Then the Markov chain (tx„w) with initial law Q is ergodic. 

Proof. The proof of this proposition is quite standard — the result has been 
used at various levels of explicit detail in the literature — although the general 
setting makes the use of ergodicity of P a bit subtle. Kozlov [89] proves this by 
way of a functional theoretical argument; we will follow a probabilistic argument 
from Berger and Biskup [11]. 

Let A be the event on the space of trajectories (wn)„gz that is shift invariant. 
Explicitly, if 9 is the Markov shift, (6'aj)„ — ujn+i, we have 9~^{A) ~ A. Let 
denote the law of the trajectories induced by the Markov chain with stationary 
measure Q. Our goal is to show that /i(A) e {0, 1}. 

The first part of the proof is the classical approximation argument that drives 
the proof of more or less every known zero-one law. Define the function 

fiiuo) E^{lA\iOo). (2.12) 

We claim that = / /i-a.s. To this end approximate A by a sequence of events 
A„ G (t(w_,i, . . . ,u!n) SO that 

WU-iAJmp) 0. (2.13) 

The shift invariance of A implies that the same holds for A„ replaced by 6" (An) 
and by 0-"(A„). 
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Invoking the general fact E ^\E ^{g\uj[))\ < Efj_\g\ and applying (2.13), we thus 
have 

\\f-E,{le..M\\,.,^, „^0- (2-14) 

Similarly, replacing 1a by Ia^a and approximating the first indicator by lgn(^An) 
and the second by Ig-n^^^-) we obtain 

\\f - E,,{^g-^^A„)hHAJ^o)\\LH^.) 0. (2.15) 

But 6'"(yl„) G ct(wo, • ■ • 7 W2ti) and 0~"(A„) G CT(w_2n, . . . , cjq) have only one 
coordinate in common and so, conditional on ujq, they are independent. This 
means 

Efj.ile-"{A„)'^e"{A„)\i^o) = -^A'(le-"(A„)|wo)£'p(l6/"(A„)l'^o)- (2-16) 

Passing to n — > oo, the right-hand side tends to /(wq)^ in L^{n) thus proving 
that f ^ P /i-a.s. 

The second step is more subtle. Indeed, we claim that /(t^uj) = f{uj) for 
all X € Z'' and Q-almost every uj. To this end let us note that, by the 6- 
invariance of A, if loq is the initial configuration of a path in A, then also uji is 
the initial step of a path in A — namely, the shifted path! A moment's thought 
shows that this implies /(i^o) ~ fi^i) M"a-S. and thus 

f{TxM = /(^) (2.17) 

for Q-a.e. ui and P°-a.e. trajectory (X„) of the Markov chain. The conditions 
on P guarantee that for P-a.e. uj, with positive probability (X„) visits any given x 
and so we must have /(txUj) ~ f{uj). The event {/ ~ 1} is thus shift invariant 
and so P(/ = 1) e {0, 1}, by the ergodicity of P. Then 

^l{A) = E^if) ^ Qif = 1) e {0, 1}, (2.18) 
where we used that the wo-marginal of is Q and that Q P. □ 

2.2. Vanishing speed 

The conclusion of Lemma 2.1 and Proposition 2.3 can be formalized in multiple 
ways. E.g., we thus know that for any / = /(w) with Eq\f{io)\ < oo, 

^ n— 1 

lim -^/(rx,c.) = i?Q/(a.) (2.19) 

A:=0 

for P-a.e. uj and P^-a.e. path (X„). But since the convergence comes from the 
Markov shift, we are not limited to functions of only one argument. Thus, 
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for instance, we also know that for any function / — f(uJo-,uJi) such that 
E(iE°\f{uj,Tx,uj)\ < oo, 

^ 71 — 1 

hm - firx.^.rx,^.^) = E^E°^f{u;,rx,u;) (2.20) 



n— >oo Tl 

k=Q 



for P-a.c. uj and P^-a.e. path (X„). This permits us to prove: 

Theorem 2.4 (Vanishing speed). Suppose P obeys assumptions (1-3) in Pro- 
position 2.3 and 

¥.(^LOo,Ax\^ < (2.21) 
Then for P-a.e. u and P^-a.e. trajectory (Xn), 

lini — = 0. (2.22) 

n— >oo n 

Proof. Our key problem is to represent as an additive functional of the 
Markov chain (t^^w). This can be done easily under the assumption that the 
environment is not periodic: 

P{{uj: T^,oj ^oj}) =0, .TT^O. (2.23) 

(Clearly, if the environment is periodic in some direction, there is no way for the 
walk to "notice" its motion through it when it makes a step in that direction.) 
We will thus prove the theorem only in this case leaving the periodic cases — 
which for ergodic P are a.s. events — to a (simple) Exercise afterwards. 
We claim that, under (2.23), we get 

n-l 

X„ = YfiTx,uj,Tx,^,uj) where /(w, w') := ^ zl{,^'=r,<..}- (2.24) 

fe=0 2 

Indeed, for almost every environment and any path of the chain, at most one of 
the indicators in the definition of / will be non-zero, and it is precisely the one 
that relates uj' to the shifted configuration oj. 

We will now apply the conclusion (2.20), but to get the conclusion of the 
theorem we need to show that 

SQS°|/(u;,rxi^)| <oo and EQE"^fiLj,Tx,Lo)^0. (2.25) 

This is a matter of a straightforward calculation. First, 

EQE°\f{uj,Tx,oj)\=EQE°\X,\ = 1e(^^l.o..|.t|) < oo. (2.26) 
Second, the absolute summability we just showed permits us to write 
E^E^fiuj^Tx.io) = EqE°X, = ^e(J2'^o,.x) 

= ^'Y'&{UJQ^:,X + UJQ^^^{-X)). (2.27) 
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To see that the last expectation vanishes, recall (2.6) to see that Ewo,^; = 
E(a;o,-:.). ' □ 

The minor trouble with periodic configurations disappears if we encode the 
sequence of environments along with the corresponding (next) step of the walk. 
This is an approach that was taken in Kozlov [89]; however, the above works 
just as well. Indeed, we pose: 

Exercise 2.5. Consider a product law on configurations iuj^a) where lo is sam- 
pled from P and a = {cFx)xeZ'' ore i.i.d. (non- degenerate) Bernoulli. Show that 
this law is ergodic with respect to the Markov shift 

(a;,(T) !-> {TXiUJ,TXiCr) with Xi sampled from P°. (2.28) 

Find a function of the (joint) environment which encodes Xn as an additive 
function of two consecutive environments. Use this to conclude that (2.22) still 
holds for almost every path of the Markov chain over u, regardless of whether 
the aperiodicity condition (2.23) holds or not. 

Notice that, for a shift-invariant configuration tu, the condition (2.21) reduces 
exactly to the first moment condition in the SLLN. So (2.21) should generally 
fail once (2.21) is violated, although exact conditions under which this is true do 
not seem to be available. The same should apply (under a different condition) 
when only convergence in measure is in question. 

The following lemma, which arose in the writing of a proof in Biskup, Louidor, 
Rozinov and Vandenbcrg-Rodes [16], can sometimes be useful in applications: 

Lemma 2.6. Let f G LlogL(Q) and suppose that P obeys assumptions (1-3) 
in Proposition 2.3. Then for P-a.e. a;, 

hm -EiCy. /(^^.^)) = ^q/(^)- (2-29) 
In particular, the limit exists P-a.s. 

Proof. Without loss of generality assume that / > and recall that ilogi(Q) 
is the space of functions / such that /log|/| S L^((Q). By Wiener's Domi- 
nated Ergodic Theorem (e.g., Petersen [116, Theorem 1.16]) these functions are 
distinguished by the fact that 

^ n— 1 

r :=sup- ^/orx, eL^(Q). (2.30) 

Since i X)fe=o f°'''Xk ^^'C dominated by f* and tend to £'q/(w) P°-almost surely 
for P-a.e. ui, the result follows by the Dominated Convergence Theorem. □ 

A subtlety of the above statement is that although the averages ^ X]fc=o / ° 
converge almost surely and in L^{Q^ P°), this is not enough to guarantee 
convergence in L^(P°), for P-a.e. w. A useful step towards understanding this 
is solving: 
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Exercise 2.7. Construct a sequence of random variables Z„ such that Z„ — > 
Z almost surely in , but such that, for some a -algebra , the conditional 
expectations E(Zn\£/) do not converge almost surely. 

The above arguments are useful even for the continuous-time versions of our 
random walk. Indeed, we can combine Exercise 1.10 with Theorem 2.4 to solve: 

Exercise 2.8. Let P be ergodic with P(7rj^(0) > 0) = 1 and ^tt^{Q) < oo. Then 
for¥-a.e. lo, the VSRW does not escape to infinity — i.e., no blow-ups occur — 
in finite time. 

2.3. Martingale (Functional) CLT 

Once a variant of the Law of Large Numbers has been established the next natu- 
ral question is that of fluctuations. In order to discuss all aspects of this question 
in a reasonably pedagogical fashion, for a while we will restrict attention to a 
class of toy models in which the environment has the following properties: 

Assumptions 2.9 (Toy- model assumptions). For some a G (0, 1) andP-almost 
every lo, 

(1) LOxy = unless \x ^ y\ ~ 1 (nearest-neighbor environment). 

(2) For each coordinate vector and each x £ 11^ , '^x-ci,x = ^x,x+ci- 

(3) a < ujx,x+ci < ^ for all i and all x. 

In other words, the environments are nearest-neighbor, elliptic and the conduc- 
tances are constant along the edges on each line of sites in . 

What makes these environments special is: 

Lemma 2.10. Let := cr(^oi ■ ■ ■ ■, ^n)- For all environments above, {Xn, -^n} 
is a martingale. 

Proof. Any environment satisfying conditions (1-2) above has the property that 
the local drift, 

]/(w) := E^JX,), (2.31) 

identically vanishes. To see how this implies the claim we note that, by the 
Markov property the law of Xn+i — A„ conditional on A„ is that of Xi in 
distribution P^^ ^. Hence, 

£;S(A„+i|^„) = A„ + ViTx^iu) = A„ (2.32) 

and so A„ is a martingale. □ 

We remark that more general (particularly, non-reversible) cases of such bal- 
anced environments have been treated by Lawler [92], Guo and Zeitouni [77] 
and, quite recently, Berger and Deuschel [13]. The main issue dealt with in 
those papers is a construction, and proper control, of an ergodic, invariant law 
on environments. 
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Returning to the setting of Toy Models, the fact that X„ is a martingale 
with bounded increments immediately implies, via Azuma's inequality, Gaussian 
bounds on its tails. Explicitly, for any unit vector e S M'' we will have 

P°(e-X„ > AV^) <c-^'/^ (2.33) 

However, to get the desired CLT we will have to invoke a more delicate tool 
which is: 

Theorem 2.11 (Martingale Functional CLT). Let {M„,^„} he an R-valued, 
square-integrable martingale such that the following conditions hold: 

(LFl) There is € [0, oo) such that for all t > 0, 

[tn] 

lim -Y E{\Mk+i-Mknn) ta^ (2.34) 

in probability. 
(LF2) For each e > 0, 

1 " 

lim -^i?(|Mfc+i-A4|2l|^, — > (2.35) 

fc=0 

in probability. 
Then for each T > 0, the law of 

t ^ ^(A^Lt^J + - [tn\){M^tn]+i - A^L*"] )) (2-36) 

on C([0,T], tends to the Wiener measure with EBt = and EBf = tcr^ 

This is what is sometimes referred to as the "Lindebcrg-FcUer Functional 
CLT," although this is only thanks to the formulation which is borrowed from 
the context of sums of independent random variable (the Lindeberg-Feller CLT, 
see, e.g., Durrett [51]). The result for martingales is, in this formulation, first 
due to Brown [23]. Derriennic [46] gave a thoughtful survey of these results; 
unfortunately, the full version of his paper is somewhat hard to get hold of. 

A simple way how to understand the scaling of the martingale paths to Brow- 
nian motion is via Skorohod embedding. Explicitly, we have: 

Theorem 2.12 (Skorohod [125], Strassen [127] and Dubins [49]). Suppose that 
{Mn,.^n} is a square-integrable (real-valued) martingale with E{Mq) = 0. Then 
there is a sequence of integrable stopping times (Ti) with Tq — and Ti+i > Ti, 
such that 

Lawof(M„)„>o = Law of (BT„)n>o- (2.37) 

The history of this result is roughly as follows: Skorohod [125] noted its va- 
lidity for sums of independent random variables, Strassen [127] observed that 
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Fig 2.1. An example of the Random Conductance Model satisfying Assumptions 2.9(1,2), 
but not the uniform ellipticity requirement in part (3). Here, for each horizontal or vertical 
line of edges in T? , we independently retain, resp., drop all edges with probability p, resp., 
1—p. The resulting random subgraph ofT? is almost-surely connected and the conclusion of 
Lemma 2.10 holds for almost every sample thereof. 

it holds even for martingales and Dubins [49] finessed an important technical 
detail where the construction of the stopping times can be done purely on the 
path-space of the Brownian motion (i.e., without reliance on additional ran- 
dom variables). 

Returning to the above Martingale CLT, condition (LFl) guarantees that 
Tn/n -> which means that the time change between the martingale and the 
Brownian motion is asymptotically linear. The condition (LF2) ensures tightness 
in the space of continuous paths (i.e., the Brownian motion will not wiggle 
too far from the piece-wise linear path interpolating the martingale values). 
The Skorohod representation only applies to R-valued martingales, hence our 
restriction to those. 

We can now finish the proof of: 

Proposition 2.13. For any shift- erg odic environment law P satisfying (Toy 
Model) Assumptions 2.9 and forP-a.e. sample from it, the law oft i— > 
linearly interpolated into a continuous path, tends to Brownian motion. 

Proof. We already know that Xn is a martingale for the filtration ^„ :— 
(t{Xo, . . . ,Xn) so we need to verify the conditions of the above theorem. This 
will be done again by using the point of view of the particle. By the Cramer- 
Wold device it suffices to prove the convergence for the projections onto all 
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vectors in W^. Fix a vector e G and consider the function 

f{u):=El{\e-X,\'') (2.38) 

and define M„ := e • X„. The Markov property guarantees 

El{\Mk+i - Mu?\.^k) = firx^io) (2.39) 

and since / is bounded and the environment is elliptic, (LFl) follows with 
(7^ Eqf{ijj) by (2.19). The condition (LF2) is trivially satisfied and so we 
have the result. □ 

Notice the (somewhat counterintuitive) fact that we prove a CLT type of 
result by invoking a LLN type of result. But this is not so strange when we realize 
that for convergence to Brownian motion we need three things: asymptotically 
independent increments, their zero mean/second-moment property and their 
stationarity. The former two properties can be safely attributed to the use of 
martingales, but for the last one — and, in this setting, the most difficult one 
— we need to use the Ergodic Theorem and thus the machinery originally 
developed for the LLN. 

Exercise 2.14. Consider the example of a random environment in Fig. 2.1. 
Show that, for almost every realization of this environment, the Martingale CLT 
applies. Characterize the variance of the limiting Brownian motion. 



2.4- Martingale approximations and other tricks 

The derivations in the preceding sections, however elegant, hinge on the crucial 
assumption of vanishing drift. Unfortunately, this is not what one can (and 
wants) to ask from a generic Random Conductance Model. Historically, this 
puts us somewhere in the first half of 1980s when people made first successful 
attempts to address the CLT in this level of generality. We will follow Kipnis 
and Varadhan [86] where the following strategy was taken: 

(1) Represent X„ as the sum of a martingale and an additive functional of (a 
single state of) the Markov chain on environments. 

(2) Approximate the additive functional by a martingale with an error that 
can be controlled at the level of the CLT. 

The first step can be achieved trivially: 

X„ = [Xk+i -Xk- E(Xu+i - Xu\,^k)\ + ^i^k+i - Xkl^k) (2.40) 

fe=0 fe=0 

The first sum on the right is clearly a martingale — call it Af„ — while E{Xk+i — 
Xk\^k) = VijXk^) makes the second part an additive functional of the Markov 
chain (TXf.ijj). (Note that we already know that Xn as additive functional of two 
consecutive environments, but for the application of the Martingale Functional 
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CLT the dependence on a single environment is much easier.) Now we need 
to write 

n-l 

Y,V{tx,u)^M;^ + E^, (2.41) 

fc=0 

where maxfc<„ \Ek\/y/ri tends to zero in probabihty. This can be done under 
proper conditions but one then faces the (rather extreme) difficulty that M„ 
and Afj'j are not independent. 

To sec how an additive functional of a Markov chain can be approximated 
by a martingale, consider a Markov chain on a state space fl with transition 
kernel Suppose g: 57 R is a function such that g G Ran(id— In other 
words, we require 

g = h~ (2.42) 

for some function /i : O — > M. If wq, wi, . . . denote the successive states of the 
Markov chain, then a similar trick to the one used above yields 

n—1 n— 1 

^ g{uk) = h{uo) - h{uj,,) + ^ [h{uk+i) ~ ^h{ujk)] ■ (2.43) 

Set En := h{ojQ) — h{ujn) and define Af,j to be the sum. By the Markov property, 

^h{uk)^E{h{uk+,)\a{uo,...,ujk)), (2.44) 

which implies that (M^) is a martingale. Of course, in order to have a useful 
statement, we need that this martingale is properly integrable, which means 
that the Poisson equation (2.42) must be solved with h in, say, L^. As we will 
comment in a minute, this may be quite a challenge to prove (and in fact, it is 
often too much to ask). However, such considerations are entirely unnecessary 
for finite-state Markov chains: 

Exercise 2.15. Consider a Markov chain with a finite state space 17 and a 
stationary measure Q. Let 5: 17 — > K. satisfy Eqg = 0. Show that, for Q-a.e. 
initial state wq, the law of 

tends to a mean-zero normal random variable. Characterize its variance. 

This statement is actually one of the main results of a note due to Gordin 
and Lifsic [69] . It will be easy to see that the result generalizes to arbitrary state 
spaces under the condition that g € Ran(id — ^) — which we take to mean 
that (2.42) has a solution h £ L'^{Q); the error, 

E„:=h{ojo)-h{ujn), (2.46) 

is then trivially bounded in L^. However, a bounded error is a luxury that we 
do not need; indeed, for the purpose of the CLT one can tolerate errors up to 
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o{y/n) — particularly, if that brings the benefit of weaker conditions on g. A 
milestone achievement in this vain is the result of Kipnis and Varadhan [86] who 
proved the following theorem: 

Theorem 2.16. Suppose that a Markov chain on state space f2 with transi- 
tion kernel 0^ is reversible with respect to Q. Consider the law on trajectories 
{^n)n>t) where ujq is sampled from Q. Let g € L'^{Q) with Eqi^g) = 0. Then the 
law of the (2.45) tends to a (zero-mean, finite-variance) normal random variable 
A/'(0,(Tg) if and only if g £ Ran([id— ^]~^/^) or, equivalently, 

sup Eq(g{u){l + e - ^y^g{uj)) < oo. (2.47) 

Moreover, the supremum equals and the convergence extends (with the limit 
given by Brownian motion) even to paths (linearly) interpolating the values of 

t^n-y^Y^tigi^u). 

Note that the claim concerns the averaged law; no statement about a typical 
starting point is made. This is one of the deficiencies we will have to address in 
detail when proving the quenched invariance principle in the next two sections. 
The original method of proof in [86] was to consider the spectral measure 
associated with the function g and the operator 0^ on L^(Q). This measure has 
the property that, for any F G L^(^g), 

(5,i^(^)5)L^(Q) = / ^(A)A^s(dA). (2.48) 

The Kipnis- Varadhan condition (2.47) can then be written as 

^^(dA) < oo. (2.49) 

Notice that the spectrum of 0^, and thus the support of fig, is contained in 
[-1,1]- 

Exercise 2.17. Show that if (tUn) is a stationary Markov chain with loq dis- 
tributed according to Q, and g G L^(Q), then the variance of (2.45) tends to the 
quantity in (2.49). 

The spectral measure is a very interesting object in its own right due to 
the connection with the area of random Schrodinger operators. What is quite 
puzzling is that we do not have any substantive information to report on: 

Problem 2.18. Describe the connection between the spectral properties of the 
generator of the random walk among conductances to — many of which, as 
is well known, are same for a.e. uj — and the generator ^ := 0^ — id of the 
Markov chain on environments. 

Let us make some remarks on how the history of the above ideas seems 
to have evolved. First, the idea to decompose additive functionals (of general 
stationary ergodic processes) into a martingale and an error is presumably due 
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to Gordin [68] who also had the insight to characterize the objects in terms of 
their functional-analytic (rather than mixing) properties. Gordin and Lifsic [69] 
then applied this idea in the specific context of finite-state Markov chains. 

The understanding that martingale approximations can be the ultimate pas- 
sage to limit laws for random walks in random environment seems to have 
grown out of the work of Papanicolaou and Varadhan [112]; the predecessors 
of this work were mostly focused on periodic environments. An alternative ap- 
proach based on resolvent methods was devised by Kiinnermann [91]. The above 
(Kipnis- Varadhan) Theorem 2.16 more or less closed the matter for the annealed 
law in reversible cases. Two natural ways to generalize Theorem 2.16 arc as fol- 
lows: One is to go beyond the annealed law and the other is to extend beyond 
reversible Markov chains. Both of these directions are far from settled and both 
constitute a subject of intense research. 

We will expound on how to go from annealed to quenched laws in the rest 
of these notes. Concerning departures from reversible situations, two lines of 
thought are generally being followed: One approach, drawing on the functional- 
analytic ideas, goes by imposing (and checking) various sector conditions (e.g., 
OUa [110], Sethuraman, Varadhan and Yau [121], Horvath, Toth and Veto [80]). 
The role of these conditions is to control the antisymmetric ( "non-reversible" ) 
part of the generator by the symmetric one. Another approach goes by imposing 
decay-rate conditions on time-correlations (e.g., Maxwell and Woodroofe [103], 
Derriennic and Lin [47], Peligrad and Utev [113], Klicnarova and Volny [88], 
Volny [133], etc.). However, unlike the reversible situations, it does not seem 
likely that a single condition will eventually cover all cases of interest. 

3. Harmonic embedding and the corrector 

Although the subject of martingale approximations is very attractive and use- 
ful, in the sequel we will adopt a different approach that emphasizes the geo- 
metrical component of the problem over its analytical component. To motivate 
this approach, consider the explicit example of the simple random walk on the 
two-dimensional supercritical percolation cluster. When the local drift V{u!) is 
non-zero, then this is because there is an odd number of neighbors of the origin 
and the origin thus no longer lies in the bary center of its neighbors. The mar- 
tingale defect can therefore be thought to arise from the use of the geometric 
embedding of the graph, before the edges got removed. 

This suggests an idea that one might instead try to look for a different, 
harmonic embedding for which V would trivially vanish. A moment's thought 
shows that such an embedding is easy to find in any finite box using a com- 
puter — just freeze the positions on the boundary and then ask the computer 
to sequentially pass through all vertices and always put them at the center 
of mass of their (graph-theoretic) neighbors. It turns out that this procedure 
rapidly converges and leads to a picture as in Fig. 3.1. How such an embedding 
is generated without recourse to finite volume is a slightly more complicated, 
although not unsolvable problem. The main new ingredient will be the reliance 
on homogenization theory. 
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Fig 3.1. Left: A sample of a percolation cluster at 50 x 50 square with bond probability 
p = 0.65. Only vertices connected to the boundary are retained. Right: The corresponding 
harmonic deformation obtained by relaxing the positions (except those on the boundary) to 
make each vertex lie in the "center of mass" of its (graph-theoretic) neighbors. Note that all 
dangling ends — parts of the cluster attached only by one edge — collapse to a point while 
the components attached by exactly two edges line up along a linear segment. 



Here and henceforth we wih make repeated use of this notion: 

Definition 3.1. We will henceforth say that P obeys the "usual conditions" if 
it satisfies the conditions (1-3) in Proposition 2.3. 

These are exactly the conditions that guarantee the existence and ergodicity 
of the Markov chain on the space of environments. 

3.1. Minimizing Dirichlet energy 

We begin with some motivational observations for general reversible Markov 
chains that will explain in more detail how Fig. 3.1 was generated. Suppose a 
countable set Y is given along with the collection of (non-negative) conductances 
i'^xy)x,yi£y subjcct to restrictions (1.1-1.2). Suppose in addition the irreducibil- 
ity condition: for each x,y G y, there is an n > with P^(x, y) > 0. For a finite 
set ^ C we then define 

{x,y}nA^(/) 

to be the Dirichlet energy in A for the potential /. The following is well known: 

Lemma 3.2 (Dirichlet principle). Let A C be a finite set with Y \ A ^ % 
and let g : f ^ M be a bounded function. Then the infimuni 



inf{£:^(/): fy^A = 



(3.2) 
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is achieved by the unique solution to the Dirichlet problem 

('U/(x)=0, ifxeA, 
\f[x)=g{x), zfxer\A. 

Proof. Pick X ^ y and any function /. Let fx be defined by 

/.W-j^^^!' (3.4) 

We claim that whenever x G A, the "move" f ^ fx demonstrably lowers the 
Dirichlet energy, £A{f) > ^Aifx)- This is seen from the identity 



y 



^xv[f{y) - fix)]'^ 

^J2uJxy[fiy) - P./(x)]' - P./(x)]', (3.5) 

y 

which is proved by optimizing the left-hand side over possible f{x) — this shows 
that the minimum is achieved at Pi^f{x) — and using that, for h{x) := Ax^ + 
Bx + C with A > 0, if .Tniin is the minimizer then h{x) = h{xmin) + A{x — Xmui)'^ ■ 
The explicit control of SaIJ) — ^Aifx) shows that, applying the averaging 
f{x) !->■ Pu:f{x) keeps lowering the Dirichlet energy as long as f{x) ^ Pw(/). 
Furthermore, 

inf f{z) < inf fx{z) < sup fx{z) < sup f{z), (3.6) 

zer zsr ^e-r z^~v 

and so any minimizing sequence of £A{f) in (3.3) is bounded. Reducing to 
subsequences if needed, we extract a limit which then obeys / = P^jf on A and 
thus solves the Dirichlet problem. To see that the solution is unique, note that 
f = P^f on A implies that / cannot have (strict) local extrema inside A. In 
particular, we have the maximum principle: 

inf (^(z) < min/(z) < max/(z) < sup giz). (3-7) 

Linearity guarantees that the difference between two solutions to (3.3) solves 
(3.3) with g := 0. The maximum principle ensures that the difference must be 
zero. □ 

The above proof suggests that we could perhaps use the Dirichlet energy as 
a kind of measure of distance from a harmonic function. We will explore this 
very soon in a more general context. However, the argument also highlights a 
difficulty associated with attempts to "harmonize" the linear function f{x) = x 
in infinite volume. Indeed, the full-lattice Dirichlet energy of such an / is infinity 
and so the procedure does not make sense. 

This problem is not unknown from other situations and it naturally leads us 
to a guiding principle of homogenization theory: Instead of trying to find the 
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Fig 3.2. A graphical solution of the electrostatic problem depicted in Fig. 1.3 based on a 
harmonic deformation of the underlying graph. The electrostatic potential changes linearly in 
the height (more precisely, the y-coordinate) of the point. In particular, the potential at the 
vertex marked by the star — originally, the origin of coordinates — is proportional to the 
ratio between its distance to the top and the bottom plates. 

deformation of the linear function f{x) = x that is harmonic with respect to 
at all locations for one given a;, we will solve the problem at one specific 
location — namely the origin — but simultaneously for all w. Technically, this 
amounts to replacing the space £'^{iTi^) associated with the Markov chain (X„) 
by the space L^{Q) associated with the chain (rx^w). The advantage of working 
on L^(Q) is that, unlike tt^:, the measure Q is finite. 

3.2. Weyl decomposition and the corrector 

To motivate the forthcoming definitions, recall that the process of substituting 
^uf{x) for /(x) applied to the function f(x) := x would replace the value x 
hy X + ^ [Xi ) . From the point of view of the particle it makes sense to shift 
this so that the origin of coordinates will not be moved under this action and 
so we may in fact want to replace xhy x-\- ^ ^S(^i)- The difference 

E^ cj(^i) ~ ^S(^i) is in the form of a gradient, 



We are thus led to minimizing the functional 



(3.8) 




(3.9) 



X 
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over all, say, local functions tp — (p{uj). Here we recall that ip — (p{uj) is said to 
be local if it is a bounded, continuous function of a finite number of luq^x^s. To 
see how homogenization translates finite- volume quantities to functionals over 
space of environments, it is instructive to solve: 

Exercise 3.3. Consider the Dirichlet energy £ajv(/) from (3.1) for the set 

Atv := [— iV, A^]'* n Z**. Fix a local function ip ~ fi^^) and set f to x + Wip 
defined by {x + Wip){z) :^ z + (p o Tz{uj) for z G Z"^. Assuming u is a sample 
from an ergodic measure P, carefully check that 

lim -^£A^{x + Vip) =E(y2ujo,,\x + \7M^)\]- (3-10) 

N^ac |Ajv| / 

For technical reasons it will be advantageous to interpret (3.9) as a quadratic 
form on vector fields. Let J\f denote the set of admissible jumps of the Markov 
chain, 

Af := {x e Z'^ : V{uja,^ > 0) > 0} U {0} (3.11) 

By a vector field we will then mean a (measurable) map u: x A/" — > M.'^, i.e., 
a vector valued function u — uicj, x) indexed by environments and points in Af. 
We will always set 

u(a;,0):=0 (3.12) 

by definition. 

An example of a vector field is a potential field Vt/? where \/ip{uj,x) := 
Va;</5(w). Any potential field is curl-free in the sense that it obeys the cycle 
conditions. These conditions state that for any sequence ccq, xi, . . . , x„ := xq of 
vertices in Z'^ such that x^+i ~ Xi E Af for all i we have 

n-l 

^u{Tx.u,Xi+i - Xi) = 0. (3.13) 

1=0 

In light of our convention (3.12), whenever J\f generates all of Z'' (as an additive 
group), this turns out to be equivalent to 

u(u), X + z) — u(uj, x) ~ u^TxUj, z), whenever x, z,x + z €z Af, (3.14) 

The vector fields that obey this property (for all cj) will be called shift covariant 
(sometimes they are called stationary). Note that from (3.14) we automatically 
have u(lo,x) = —u{TxitJ,-'x). 

As already alluded to, all potential fields are shift covariant. Another example 
of a shift-covariant field is the position field, x{u},z) := z. As we shall see later, 
the position field and the potential fields generate the vector space of all shift- 
covariant fields. The reason for singling out shift-covariant fields is that they 
correspond to gradients of lattice functions. The following exercise details this 
connection: 
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Exercise 3.4. Assume the irreducibility condition P(sup„>]^ P"{0,x) > 0) = 1, 
for all X G IfK Show that for any shift covariant u there is a (F-a.s.) unique 
function U ~ U{uj,x) with x e Z'' such that 

J7(w,0) = (3.15) 

and 

) - U{uj, x) = u{t^uj, z), X eZ"^, ze TV. (3.16) 

To indicate that the vector field u(u;,x) and the function U{uj,x) are related 
as in (3.15-3.16), we will sometimes write u = gradt/ or say that U is an 
extension of u to Z''. 

The expression 

{v,w) ■='E{ ^ wo,xw(w,x) • «;(i:j,x) j (3.17) 

defines a natural inner product on the set of vector fields; the dot in w(aj, ■ 
t«(a;,x) stands for the usual (Euclidean) dot product in . This inner product 
defines a natural L^-norm; a minor technical problem — which has often been 
overlooked in the literature — is that (w, u) = does not imply that m = 0, 
only that u!q xu{uj,x) = for all x € M. A standard approach would be to 
factor the space of vector fields by the equivalence relation u u' whenever 
{u — u',u — u') = 0. However, this is unnecessary once we restrict attention to 
shift-covariant fields (and impose a proper non-degeneracy condition). Indeed, 
define the set 

^cov shift covariant, {u,u) < oo} (3.18) 

and set ||u|j/,2^^ := {u,u)^^^. It is then not too hard to solve: 

Exercise 3.5. Assume the irreducibility condition P(sup„>2^ P"(0,x) > 0) = 1, 
X E Z''. // 11^11^2 = then u{u!, x) = for all x G J\f and P-a.e. uj. 

Once the L^-structure is in place, we note that potential fields define a natural 
closed subspace 

Ll ~ {V^: if local} (3.19) 
of icov With this space comes the orthogonal decomposition 

Ll^, = Ll(B{Ll)^. (3.20) 

It turns out that the vector fields from (iy)^ '^''^^ quite well characterized. 
To see that explicitly, define the divergence div(wM) by the formula; 

div(cJM) := ujQ^x [u{u!, x) — uir^uj, —x)] (3-21) 

X 

where the bracket simplifies to 2?i(aj,x) once u is shift covariant. Thinking of 
u{ui,x) as the flux from to cc, the first term on the right (including wq.x) 
corresponds to the total flux out of the origin and the second one to the flux 
into the origin. 
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Lemma 3.6. Foru^ -^cov; we have u £ (-^v)^ if and only if div (uju) —OforP- 
a.e. Lu. In particular, ifU is a function such that u = gradt/, then L^j U{lu,x) = 
at all X and P-a.e. w. 

Proof. Pick a local function ip and note that 

u(T_^i:j,a;) - wckmIwjX)]^ (3.22) 

= — E((y9(aj)div(ti;it)), 

where we used u G L^^v split the second expectation into two terms and then 
relabeled x for —x. It follows that if {u, \/(p) ~ for all local functions, then 
div(a;u) = P-a.s. and vice versa. 

For u = grade/, a simple calculation shows div(cjw) = 2\-^ U. With the help 
of shift covariance the condition div(wit) = then forces L^^ U{uj, •) = 0. □ 

Lemma 3.6 shows that the fields in (Ly)"*" are, after multiplication by oj, 
necessarily divergence-free — and are thus sometimes referred to as solenoidal 
fields. The orthogonal decomposition (3.20) is thus an analogue of the Weyl 
decomposition from differential geometry. For readers familiar with basic elec- 
trostatics, the function U — associated to a shift-covariant field u — can be 
thought of as an electrostatic potential while uju plays the role of an electric 
current. The fact that potential difference and current are related by way of a 
multiplication by w is a demonstration of Ohm's law of electrostatics. See Doyle 
and Sneh [48] and/or Sect. 6. 

A natural next question to ask now is whether there are any solenoidal fields 
at all. For nearest-neighbor, constant conductances, a perfect candidate for a 
solenoidal field is the position field which simply assigns x{lu,x) := x. (Indeed, 
this function is discrete harmonic with respect to the homogeneous Laplacian 
on Z'^ and so it obeys the conclusion of previous lemma.) Of course, once the 
conductances are not constant, div(cjx) — which equals twice the local drift 
V{Ld) — is generally non-zero, but one can still hope that x has a non-trivial 
projection into the subspace (Ly)-'-. This is all expressed in: 

Proposition 3.7. Suppose P obeys the "usual conditions" and, in addition, 
assume that 

Ef ^a;o,:r|a;n < oo. (3.23) 

Then there is a function ^ = '${uj,x) defined for all x £ Z'^ with the properties: 

(1) Harmonicity: L^^'^ [uj , x) = for all x G Z'' and P-a.e. uj. 

(2) Shift covariance: 5'(ti;,0) = and 



*(cj, x + z) - "^{uj, x) = ^{t^lo, z), X, z € Z''. 



(3.24) 
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(3) Square integrability: EqE^\^{uj,Xi)\'^ < oo. 

In addition, for any minimizing sequence ipn of the function (3.9), we have 
V(y5,i — > x(w, a;) in L^^^ where x the corrector that is given by 

X{uj,x) ■.= 'f{u,x)-x. (3.25) 

Furthermore, ^! = 5'(a;,x) with x restricted to Af is the orthogonal projection 

*(w, •) := proj(i2^)_Lx(a;, •)■ (3.26) 

The infimum of (3.9) over all 'Vip e Ly is exactly ||5'||^2 • 

Proof. The proof could be simply started by defining 5* via (3.26) and then 
checking the stated properties based on facts from the theory of abstract Hilbert 
spaces. However, it will be more instructive to prove some of the those claims 
directly in the present setting. 

First note that the object in (3.9) can be interpreted as 

e(^^cjo,.|x + v,</j(cj)|'^ = ll^ + ^'^lli?„v- (3-2^) 

The condition (3.23) then guarantees that (3.9) takes a finite value for all local 
functions. Since it is also positive, we can pick a sequence (p„ for which it tends 
to its infimum. The parallelogram law then yields 



^WVifn - ^fm\\l2^^ - \\X + ^^nWh^ + \\x + V'f„,\\l2^ 



a; + V- 



. (3.28) 

The first two terms on the right both tend to the infimum while the last term is 
bounded by twice the infimum. It follows that V</3„ is Cauchy in L^^^ and so it 
converges to a vector field that we denote by x- This is the corrector in (3.25). 

Since x is a limit of gradients, it is shift-covariant and so it extends to a 
unique function on . Now we define ^' := x + x and note that ||^'|||2 is the 
infimum of (3.9). This implies that for all local functions </? and all e, 

ll^ + eVv^lli^^^ > ll^-ili^^^ (3.29) 

Expanding the left-hand side and taking e — > yields (^1', Vi^s) = for all local 
functions, i.e., 4* e [L"^)^ . By Lemma 3.6(1), ^ is Lcj-harmonic. □ 

Obviously, the conditions (1-3) in the above proposition can be satisfied by 
5" := 0; it is thanks to (3.26) that this can generally be excluded. (However, we 
could still have that 5* is identically zero; see Exercise 4.3.) A question might also 
arise whether the function is uniquely determined by the above properties. 
Biskup and Spohn [19] showed by fairly soft arguments that this is indeed the 
case. In fact, one even has a stronger statement: 

Ll^^ = L%®{A^: AeGl.{^,d)], (3.30) 
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with Ai!{uj,x) denoting the vector whose i-th Cartesian coordinate is given 
by '^jCiijej ■ where A = (fly). The position function and the poten- 

tial fields thus generate all shift-covariant square-integrable (R''- valued) vector 
fields. (Notwithstanding, see Problem 4.18 for a very non-trivial generalization 
of this question.) Quastel [117] has derived a similar result to (3.30) albeit with 
the use of Poincarc inequality and spectral-gap estimates. 

It should be emphasized at this point that the above constructions have 
been quite standard — albeit perhaps in different context and using different 
notations — in various contributions dealing with homogenization theory. An 
application of these techniques to random walk in random environment was done 
somewhat independently in the Western school by Varadhan, Papanicolaou and 
coauthors and in the Russian school by Kozlov. 

In particular, Kozlov's well-known paper [89] contains an extended version of 
the Weyl decomposition of vector fields — which he calls forms — into the sum 
of a gradient field, a harmonic field and a constant field which applies even in 
non-reversible situations. Apart from strong ellipticity, the main requirements 
for this decomposition in [89] are: 

(1) There is an measure Q which is invariant for the Markov chain on envi- 
ronments and absolutely continuous with respect to P. 

(2) The reciprocal value of the Radon-Nikodym derivative ^ is in L^(P). 

While the absolute continuity of an invariant measure is usually somewhat chal- 
lenging, it is the second condition that is invariably nearly impossible to check 
directly in any realistic (non-reversible) situation. We note that although Ko- 
zlov's paper is known to contain inconsistencies, it puts forward a number of 
good ideas and is thus a very recommended reading for anyone with interest in 
this subject. 

The construction of the harmonic deformation can be performed rather seam- 
lessly even in the case when tt^j^x) is zero at some vertices. What we need to 
assume is that there is a P-a.s. unique infinite component "^^ao of vertices with 
T^L^i^) > such that the conditional measure 

Po(-) =P(-|Oe'^oo), (3.31) 

with expectation denoted by Eq, satisfies the following variant of the "usual 
conditions" : 

(!') Po(7r^(0) > 0) = 1 (which holds trivially) and Eo7r^(0) < co. 

(2') Pq is irreducible in the sense that, for every x € Z'' with Po(x S '^oo) > 0, 

Po(w: supP^(0,x) > Olx e <^oo) = 1- (3.32) 

Tl>0 

(Condition (3) for measure P is not needed for now, the translation invariance 
of P suffices.) 
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Exercise 3.8. Suppose that ^^^o cind Pq are well defined and assume condi- 
tions (V) and (2') above. Suppose also (3.23). If Lpn is any minimizing sequence 
of the functional 




(3.33) 



show that V(p„(aj, •) still tends to some xi'^i ') -^cov ^^^^ ^'^ define \I/ = 
'^{uj^x) with X € 'looo which is harmonic with respect to L^^. 

The function ^ constructed in this Exercise is the harmonic embedding of 'rfoo 
that we discussed at the beginning of this section. A construction along the 
above Unes can be found in the paper of Mathieu and Piatnitski [101] for the 
problem of supercritical percolation cluster and in Biskup and Prescott [18] at 
the current level of generality. Berger and Biskup [11] give a construction which 
is based on the spectral representation method of Kipnis and Varadhan (see end 
of Sect. 2.4). Another way to define the corrector might be a result of: 

Exercise 3.9. Show that the limit in 

\im [EZiXn) - E'^iXn)] (3.34) 

exists and equals 5'(a;,.T) for P-a.e. u. 

It would be of much interest to find a solution to this problem without a 
recourse to the functional-analytic methods discussed above. 



3.3. Quenched Invariance Principle on deformed graph 

Let us now turn attention back to the problem of a random walk among random 
conductances. A simple consequence of the above constructions is: 

Corollary 3.10. Suppose P satisfies the "usual conditions" and, in addition, 
(3.23) holds. Define Mn ■= vE'(a;,Ar„). Then for P-a.e. lo and each T > 0, the 
law of 

t ^ -^{M^tni + [tn - [tn\){Mytr.\+i - M[tni )) (3.35) 

induced by on the space C{[0,T], tends to the Brownian motion Bt with 
EBt = and the covariance structure determined by 

-^E[{\-Btf]^E^El[[X--^{uj,Xi)f], AeR''. (3.36) 

Proof. By the Cramer- Wold device it suffices to prove the convergence in law 
for the projection of the process onto any vector. We will denote this projection 
(with some abuse of notation) also by M„ := A • '^{u),Xn). The filtration is as 
before: ^„ = ct{Xq, . . . , X„). 
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First, the L^j-harmonicity of 4* guarantees that M„ is a martingale so we just 
need to verify the conditions (LF1-LF2) of the Martingale Functional CLT. We 
will take care of both of these by considering the function 

/K(^):-i?S(|A^i|'l{|A/,|>K}). (3.37) 

Indeed, by property (3) in Proposition 3.7, /a' G L'^{Q) for all K > 0. Next, the 
shift-covariance of implies Mk+i — Mk = A . 'i'lTx^u!, Xk+i — Xk) and so, by 
the Markov property, 

E"^{\Mk+i - Mk\'liM,+,^M,\>K[n) = fxirx^cj). (3.38) 
It follows that the left-hand side of (LFl) equals 

^ n— 1 

(3.39) 

for K :— 0, while the left-hand side of the expression in (LF2) is bounded by 
this term from above as soon as n is so large that e^/n > K. 

Ergodicity of P with respect to translations ensures via (2.19) that the ex- 
pression (3.39) tends to EiQfK{uj) as n cxj. This verifies (LFl) with given 
by the right-hand side of (3.36), and it also proves (LF2) because, thanks to the 
Dominated Convergence Theorem, we have 

lim Sq/a(w) = 0. (3.40) 

AT— >oo 

The result now follows by applying Theorem 2.11. □ 

The above argument can be pushed through even in the case when the walk is 
restricted to an infinite connected component "^oa, as described above. One just 
needs to carefully check that the current proof of Proposition 2.3 still applies 
(details arc spelled out in Berger and Biskup [11]). However, later arguments 
might be seriously hampered by the fact that Pq is no longer shift invariant. 
This can be circumvented by the introduction of an induced shift. Namely, for 
each i = 1 , . . . , d, let 

0iUJ := Tn,(^u,)ii^u; (3.41) 

where 

n,{uj) := mi{n > 1: nc, S '^oo(^^)}. (3.42) 

The collection of maps {9i, . . . , 9^) defines shifts which preserve Pg and, in fact, 
make Pq ergodic. To see why these are well defined and the last property is true, 
consider the following exercise from abstract ergodic theory: 

Exercise 3.11. Let (,^,J?, /i) he a probability space and let A & ^ he such 
that ijl{A) > 0. Let t : ^ he a ^-preserving hijection and suppose that fj, is 
ergodic with respect to t. Let n^(a;) := inf{n > 1: t^^{x) € A} for each x € 
Do the following: 

(1) Show that UA < oo fj.-a.s. 
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This permits us to define 9{x) : 



(^^(x). Next: 



(2) Show that 0{A) = A fi-a.s. and that 9 preserves ) :— ^,{—\A). 

(3) Prove that fiA is ergodic with respect to 0. 

We will close this section with an exercise that illustrates the above abstract 
setting in one situation where explicit calculations are possible. 

Exercise 3.12. Suppose d = I and only nearest-neighbor conductances. Assume 
that P is ergodic with respect to the canonical shift on Z and suppose that 



defines a function satisfying properties (1-3) in Proposition 3.7. Conclude that 
the random walk (^(cj,X„))„>o satisfies the (quenched) invariance principle. 

We remark that the one-dimensional Random Conductance Model have quite 
intensely been studied, e.g., by Comets and Popov [33], Gallesco and Popov [62], 
Gallesco, Gantert, Popov and Vachovskaia [61], etc. A related problem is that of 
the random walk on random trees (with or without random conductances); see 
e.g., Lyons, Pcmantlc and Peres [98], Peres and Zcitouni [114], Gantert, Miiller, 
Popov and Vachovskaia [63]. 

4. Taming the deformation 

In this section our main goal is to finish the discussion of the essential steps 
of the proof of the quenched invariance principle. We will do this while leaving 
the most technically involved part, heat-kernel estimates, to the next section. 
Most of the material discussed here is quite standard; a possible exception is 
Theorem 4.7 which has not appeared in this generality before. 

4.1. Remaining issues 

Let us quickly review what we have accomplished so far. First, we used the 
examples of the balanced environments to isolate the martingale property as 
the key vehicle that will get us to the CLT (Section 2.3). Then, in the situations 
which are not balanced, we introduced a new embedding of — described by 
the function above — that again makes the random walk into a martingale 
(Proposition 3.7). On this embedding we succeeded in proving the convergence 
to Brownian motion (Corollary 3.10). However, two issues remained unresolved: 




(3.43) 



Verify that 




(3.44) 
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(1) The limiting Brownian motion may be degenerate to a point. 

(2) The harmonic embedding may be quite distorted from the original lattice. 

Although the answer to (1) is ultimately related to the answer to (2), we will 
first focus on (1) as it is easier. We will start by solving Exercise 3.12. 

It is easy to check that the function ^E* from (3.44) is harmonic with respect 
to L^. This follows from the calculation 

^^-(0;, x) = ujx,x+i [^{uj, a; + 1) - *(a;, x)] 

+ a;:,_i,4«'(w,a;- 1) - *(w,a;)] 



— ^x,x-\-\ r i^x-X,'!, 

The shift-covariance is a consequence of the additive form of the expressions in 
(3.44) while integrability follows from 



i?Q£;Xc.,Xi)2 = iE 



/ C \2 ( C 



1 



(4.2) 



E(wo,i)IE(lMi)' 



which is finite and positive by our assumptions. Applying the arguments in the 
proof of Corollary 3.10, ^(w,X„) satisfies an invariance principle with a non- 
degenerate limiting Brownian motion. The remainder of the Exercise is now 
embedded into: 

Proposition 4.1. Suppose P obeys the "usual assumptions" and (3.23). In 
addition, assume that 

E(l/a;o,c,) < oo, i^l,...,d. (4.3) 

Then the limiting Brownian motion in Corollary 3.10 is non-degenerate. 

Proof. We need to show that the right-hand side of (3.36) is bounded below by 
c|Ap for some c > and all A £ M''. To this end we write 



>-^^mfEf ^ uo,,[\-{x + VM^))f 

i—l ^x—±Ci 

The expectation on the extreme right now involves only edges in the i-th coor- 
dinate direction and thus effectively becomes a one-dimensional problem. To 
overcome a possible lack of separate ergodicity, let ^ be the a-algebra of 
Te^-invariant events and let Ai := A • e^. The second infimum in (4.4) is then 
bounded below by 2AfE(l/E(W(^g. j^z/^)) which by Jensen's inequality is at most 
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2Afl/^;(l/a;o,eJ- It follows that 

E„E;[[e, . *(.,X,)f] > jJ-^^ (4^5) 

By (3.23) and (4.3) we conclude that this exceeds c|Ap for some c > 0. □ 
Note that the same argument would apply whenever the set 

{x e Z'^: E(l/a;o,:r) < oo} (4.6) 

generates all of Z"^ (as an additive group). This still does not cover the case 
of supercritical percolation (which can nonetheless be covered by an alternate 
argument) so we pose: 

Problem 4.2. Give general (and natural) conditions under which the limiting 
Brownian motion is non-degenerate. In particular, if R[x, y) denotes the effec- 
tive resistivity (1.20) from x to y in environment lo, does it suffice to replace 
E(l/a;o,:,) <oo m (4.6) hyV.R{0,x)? 

Note that we have a pointwise bound i?(0,a;) < I/lo^^x with i?(0, 1) ~ I/loq i 
in d = 1 with nearest-neighbor conductances. This suggests also: 

Exercise 4.3. Suppose d ~ 1 and let V be a measure on i.i.d. positive and 
nearest-neighbor conductances such i/iai E(wo,i) < oo and E(l/aJo,i) = oo. Show 
that the infimum of (3.9) over local functions is zero. Conclude that we must 
have x(w, x) = —x. 

A proof of (an analogue of) Proposition 4.1 appeared in Kozlov [89] and in 
de Masi, Ferrari, Goldstein and Wick [43, 44]. With a bit more effort one can 
develop a variational characterization of the inverse of the limiting covariance 
matrix by minimizing a (version of) Dirichlet energy over nearly linear flows 
(Biskup [14]). This in principle allows one to numerically approximate the co- 
variance matrix with arbitrary precision from above and below. 

Approximation arguments for the diffusion constants are at the core of the 
Kipnis-Varadhan approach sketched in Sect. 2.4. Caputo and loffe [27] studied 
periodized versions of the Random Conductance Model and the convergence of 
the effective diffusion coefficient to the infinite volume object; related work in a 
continuum context can be found in Owhadi [111]. 

4-2. Sublinearity of the corrector 

Having addressed non-degeneracy of the limiting Brownian motion, we are ready 
to move to the second — and considerably more involved — issue. The important 
thing is to realize that for our purposes it would suffice to show that 

X(w, X„) = ^Sj{lu, X„) - X„ = o(X„) (4.7) 

asymptotically along a typical path of the random walk. Indeed, once we know 
that \I'(cj, Xk) — Xk = o{Xk) we can use the martingale CLT to get vE'(a;, Xk) = 
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o{^/n) for all fc < n which then implies that also Xk — 0{\/n). But then we 
will have \E'(w,Xfc) — Xk = o{^/n) for all k < n, which means that the change 
of embedding of the graph has a vanishing effect at the diffusive scale. 

A more general version of (4.7) would be to require this for all positions in 
the lattice, not just those visited by the path. In = 1, this is not hard to get: 

Exercise 4.4. Suppose that P is an ergodic law on nearest-neighbor conduc- 
tances in d = 1. Assume E(wo,i) < oo and E(l/ajo,i) < oo. Show that 

^{uj,x) - X ^ o{\x\), |.t| -> oo, (4.8) 

and prove that the corresponding random walk satisfies a quenched invariance 
principle. (Compare also with Exercise 4-3.) 

However, the situation in higher dimensions is quite more subtle. While the 
technical details of derivations in the paper of Kipnis and Varadhan follow a 
different route, their methods can be used to show: 

Theorem 4.5. Under the "usual assumptions, " (3.23) and (4.3), for each e > 0, 
£;QP°(max|x(c.,Xfc)| >eV^) 0. (4.9) 

k<n n— >oo 

This statement will imply the so called Annealed Invariance Principle, some- 
times also called a functional CLT in probability. We will choose to formulate 
this in a form of a coupling. Here we recall that, given two probability measures 
P and P', their coupling is a probability measure Q on the product space whose 
first, resp., second marginal is given by P, resp., P' . 

Corollary 4.6 (Annealed Invariance Principle). Given a path X = {Xn), a 
time t>Q and n let 

M//"^ := ^(XLt„j + {tn - [tn\){Xyt,,^+^ - Xytn\))- (4.10) 

Under the assumptions of Theorem 4-5, for P-a.e. lo, there exists a coupling Q'^ 
of the law of t ^ Wi""^ induced by P^ and a Brownian motion t i-^ Bt with 
mean zero and covariance (3.36) so that, for each T > and each e > 0, 

EQl( sup >e) 0. (4.11) 

\0<t<T / 

Proof. (Sketch) First let us note that both (4.9) and (4.11) hold equivalently 
with expectation Eiq or expectation E. (This is because Q and P are equivalent 
and the quantity under expectation is bounded.) To prove (4.11), we will use the 
fact, implied by the Skorohod embedding, that such a coupling exists between 
the Brownian motion and the analogue of 1 1— ^ Wt defined using the martingale 
M„ := ^'(w,X„). Let Z^"^^ denote the expression on the right of (3.35). Then 
we have 

EQl{ sup >e) (4.12) 

0<t<T 
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where is induced by the Skorohod embedding. As to (4.11), we note that 
sup iW^i"^ - < ^ max |x(^,^fe)| (4.13) 

0<t<T V?l fc<LTnJ+l' 

Since the event on the right does not depend on the second marginal of Qj^, we 
thus have 

EQof sup >e) <EF^"f max \x{u:,Xk)\ > (4.14) 

Vo<t<T / \k<\Tn\+l J 

which tends to zero as n — oo by Theorem 4.5. Combining (4.12-4.14) the 
resuh fohows. □ 

We remark that when the suprcmum is dropped from (4.9), we talk about an 
annealed CLT. The averaging over the invariant measure Q in Theorem 4.5 is 
not a mere technical convenience as the statement is not strong enough to infer 
(4.11) without the expectation over environment. It actually took nearly 20 years 
after Kipnis-Varadhan's result before this issue was first successfully addressed 
and a proper quenched invariance principle proved. This was done in the work 
of Sidoravicius and Sznitman [124] who realized that one can get further with 
the help of the heat kernel estimates. However, Berger and Biskup [11] were 
later able to avoid the use of these in their argument for the two-dimensional 
supercritical percolation cluster. We will present a sketch of their argument in 
a slightly more general, albeit non-percolative, setting: 

Theorem 4.7. Let d = 2 and suppose P is an ergodic law on nearest-neighbor 
conductances subject to the conditions E(wo,ei) < and E(l/wo,ei) < /'^'^ 
i = 1,2. Then 

lim — max |x(a;, a;)| = 0, P-a.s. (4-15) 

ri-i-oo n \x\<n 

Our proof of Theorem 4.7 begins by a lemma that generalizes Exercise 4.4 
— and that even in d = 1, when the conductances are no longer just nearest 
neighbor — to all dimensions: 

Lemma 4.8 (Directional sublinearity) . Suppose d > 1 and assume P is an 
ergodic law subject to the restriction (3.23). Assume E,(l/uJo,c) < oo for some e € 
{zbci : z = 1, . . . , d}. Then 

lim ^(^^^l!^ = 0, P-a.s. (4.16) 

n— f oo fi 

Before we set out to prove this, we note that there is a small technical subtlety 
that arises from the distinction between ergodicity and directional crgodicity. 
To make this distinction clearer, we invite the reader to first solve: 

Exercise 4.9. Construct a law P on nearest-neighbor conductances that is 
(jointly) ergodic with respect to translations — i.e., V{A) ~ for all A with 
Tx{A) = A for all a; € Z'' — but not separately ergodic in the sense that there 
is a set B of environments which is invariant under translations in the first 
coordinate direction and for which < f{B) < 1. 
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Proof of Lemma 4- 8. Using shift covariance we get 



n-l 



X(a;, ne) = ^ xiTk&u}, e) (4.17) 

fc=0 

We would like to use the (pointwisc) Ergodic Theorem and To-invariance of P 
to extract the limit 

/,(.) := hm (4.18) 
and prove that it vanishes P-a.s. For that we will need to establish three things: 

(1) E|x(w,e)| <(». 

(2) Ex(t^,e) =0._ 

(3) /o is translation invariant. 

The first two items will follow from the construction of the corrector. Recall 
that we are guaranteed that V(^„ — > x in L^ov — which is a kind of weighted 
L^-space. Since for (p e L°°{F), 

E\/i,(p{uj) =EipoTs -Eip ^0 (4.19) 

it suffices to show that x('iC) £ and Wsfn xi'i^) in L^. (The former 
actually follows from the latter, but we find this order more instructive.) And, 
indeed, by the Cauchy-Schwarz inequality we get 

[E|x(w,e)|]' <E(-i-)E(c^o,c|x(w,e)|') <E(^^'^\\x\\h_<^ (4.20) 

and similarly we derive 

[E|x(c.,e)-Ve¥'„(c^)|]' <e(— )||x-Vv^„||i. (4.21) 

which tends to zero as n — )■ oo because V(/3 — >■ x in L^ov 

Finally, in order to link the limit to the expectation, we also need to show 
that /c is translation invariant. To that end pick another lattice direction e' and 
note that, by translation covariance, 

x(T6'a;, nc) ^ x{^, nc) + xl^ncW, c') - x(w, e')- (4.22) 

Dividing by n, the i^-limit of the last two terms is zero and so from the above 
L^-inclusions we conclude that fc{T(,'Uj) = /c(w) for P-a.e. u. Putting all pieces 
together the claim follows. □ 

We remark that the fact that the conditions in Lemma 4.8 are the same as 
in Proposition 4.1 is not a coincidence. Indeed we have: 

Exercise 4.10 (Sublinearity implies nondegeneracy) . Show that if (4.16) holds 
for vector e, then the component of the limiting Brownian motion — constructed, 
at this point, by the Martingale Convergence Theorem — in direction of e is 
non- degenerate. 
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Fig 4.1. A figure illustrating the last part of the proof of Theorem J^.l where one needs to 
control the corrector on a component away from the good grid by the values on its boundary. 
The existence and positivity of the densities of the good vertices along coordinate directions 
ensures that the largest such component intersecting the box [— n, n]^ is o(n) in diameter. The 
argument relies on the maximum principle for x i— > '^{lu, x). 

Our next goal is to boost the directional subadditivity — which we may 
assume for both lattice directions under the conditions of Theorem 4.7 — into 
a corresponding statement over a box of side n. To this end, let us say that the 
origin is (A', e)-good in lo if for all e S {ie^ : i = 1,2} and all n > 1, 

|x(w,ne)| < K + en. (4.23) 

A point X is then called [K, e)-good in a; if is [K^ e)-good in t^lu. By Lemma 4.8 
we know that 

P(0 is (A', e)-good) — > 1. (4.24) 

A'— >oo 

It is now an exercise to show that: 

Exercise 4.11. Fix e > 0. Show that for each p E (0, 1) and for P-a.e. w there 
is a number K = K{p,Lj) < oo such that 

(1) is (AT, e)-good in u . 

(2) The density of {K,e)-good vertices on the lines {nCi'. n G Z}, i = 1,2, is 
at least p. 

These observations permit us to define a good grid as follows. Take the two 
lines {nei : n S Z}, i = 1,2, and add to them all vertices of the form niei +71.262 
with ni, 712 G Z such that either tt-iCi or 71-262 is {K, e)-good. Call the resulting 
(random) set of vertices GK,e{^^)- Then we note: 
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Lemma 4.12. We have: 

max \xiuj,x)\ <2K + 2en. (4.25) 

|x|oo<n 

Proof. Let X :~ riiCi +^202 be a vertex in Gx.£(w). This means that, e.g., niei 
is (A", e)-good in lo. Since the origin is {K, e)-good as well, we can write 



|x('^:a;)| < \xi^,x) - x(a;,nici)| + |x(niei,x)| 
< K + e\n2\ + K + e\n2\ 



(4.26) 



But |x|oo < n implies |ri,i|, |n2| < n and so the claim follows. □ 

We now know how to control the corrector at the vertices of the good grid — 
which can be made arbitrary dense — but we still have to worry about those in 
the complement thereof. An important fact is that the connected components 
of \ GK.e{^) are finite and, in fact, that any such component intersecting the 
box [— n, n]^ has diameter o(n). This can be justified by solving: 

Exercise 4.13. Consider any shift invariant, ergodic, zero-one valued process 
on TL with the densities ofO's and I's both positive. Show that in almost- every 
sample, the size of the largest consecutive block of 1 's intersecting the interval 
[— n,n] is only o{n) as n 00. 

We can now finish the proof of sublinearity of the corrector: 

Proof of Theorem J^.l. Pick x <^Tl?\ Gif_e(w) with |a;| < n. Let Cix) denote the 
component containing x. We claim that 

|x(i^,a;)| < diamC(a;) + max |x(a;,z)| (4.27) 

This is a consequence of Lt^-harmonicity of ^' and the maximum principle. In- 
deed, define the first hitting time 

T := mi{n > : X„ ^ C{x)} (4.28) 

of the complement of C{x). Then 

^{uj,x)^EZ{^{uj,Xt)), (4.29) 

which we can rewrite as 

Xico, x) = EZ {Xt -x)+EZ (x(c., Xt)) (4.30) 

But — <diamC(x) and |x(w,Xr)| < max;jga(7(j.) |x(w, z) | so the bound 
follows. 

To finish the argument, we recall that diamC(a;) = o{n) and so we may as- 
sume that n is so large that C{x) C [~2n, 2nY. In that case max^g^cj^^) |x(w, z)\ 
is bounded by the maximum from Lemma 4.12 with n replaced by 2n. We get 

\x{ijJ,x)\<2K + Aen + o{n) (4.31) 

thus proving the claim. □ 
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4.3. Above two dimensions 

The above reasoning can be boosted to cover all ergodic two-dimensional envi- 
ronments with a finite range of jumps that satisfy the condition E(l/a;o,ei) < 00. 
However, there is an inherent problem with this approach in higher dimension; 
indeed, one can still define a good grid but this grid will no longer partition 'Z'^ 
into finite components. In an attempt to adapt the argument based on (4.29- 
4.30), one thus has to worry about two things: How long does it take to hit the 
good grid and how far will Xt be from x. This can be done but (insofar) only 
with the help of heat-kernel technology. We paraphrase a theorem from Biskup 
and Prcscott [18]: 

Theorem 4.14. Fix lo such that ■nuj{x) G (0, 00) for all x and suppose x = x{^) 
is a function and 6 > is a number such that the following holds: 

(1) (Harmonicity) The function 4'(x) := x + x(x) obeys L^'^{x) — for all x. 

(2) (Sublinearity on average) For every e > 0, 

^-T ~d l{lx(.^)l>e«} = 0- (4.32) 

x\ 1 2; I < n 

(3) (Polynomial growth) 

lim max = 0. (4.33) 

Let Y = (Yt) be the variable-speed continuous-time random walk with genera- 
tor and suppose that the following estimates hold: 

(4) (Diffusive upper bounds) For a sequence 6„ = o(n^), 

Ef,\Yt~x\ 

sup max sup — r= < 00 (4.34) 

ri>l \x\<n t>b„ \lt 

and 

sup max sup f'l^PZiYt = x) < 00. (4.35) 

n>\ |a:|<n t>6„ 

Then 

lim max = 0. (4.36) 

n-)-oo |2:|<n 71 

We remark that most of the proof of this theorem goes through even when the 
variable-speed random walk is replaced by the constant-speed walk (for which 
the bounds (4.34-4.35) may be easier to prove). This is because ^(X^) is a 
martingale for both walks. The sole point where the variable speed walk seems 
to be used is formula (5.13) on page 1338 of [18]. 

In an earlier work (e.g., Bergcr and Biskup [11, Appendix A2]) the same 
conclusion as given by Theorem 4.14 could be achieved — although perhaps in 
a less transparent way — by using the full heat-kernel upper bounds of the form 

?l(x,y)<^^---\^~y\"l^. (4.37) 
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The point of reducing the heat-kernel input to the statements (4.34-4.35) is that 
these are easier to verify than the actual heat-kernel upper bounds. We also note 
that Sidoravicius and Sznitman [124] have used the heat-kernel bounds mainly 
to control the tightness of the limiting process, while here we are using it to 
control the deformations of the harmonic embedding. (Tightness follows in our 
case from the Martingale Functional CLT.) 

A key input in Theorem 4.14 is the sublinearity-on-average claim which we 
formalize as: 

Proposition 4.15 (Sublinearity on average). Suppose d > 1 and assume P 
is an ergodic law subject to the restriction (3.23). Assume E(l/a;o,c) < oo for 
all e € {±ei : i — 1, . . . ,d}. Then for each S > 0, 

i™, ^ l{lx(-.-)l>M = 0' P-a-S- (4.38) 

\x\<n 

The proof is based on the commutative structure of Z*^ and a bootstrapping 
of the one-dimensional sublinearity established in Lemma 4.8 by induction along 
dimension. Recall the notion of a good grid GK,e introduced (in d ~ 2) earlier. 
The induction argument is contained in the following deterministic "pigeon- 
hole-principle" lemma: 

Lemma 4.16. Let A„ := [—n,n]'^nZ'^, fix e > and K < oo. For each uj there 
exists a set A = A(u!) C GK,eii^) H A„ with the properties 



|A|>(2'^,y-2'' + l)|A„| where ,^ E^MllAjl 

I An I 



(4.39) 



and 

x,yeA \xiy,uj)-x{x,io)\<2d[K + e{2n+l)]. (4.40) 

Proof. We will prove this by induction on dimension. Fix cj and for G {1, . . . , d} 
define sets aI""* of the form 

d 

A';^^:=Ann Pi {x^ {xi,...,Xd): Xj =:m^} (4.41) 

for sonic m2 , . . . , <E [— n] n Z as follows: We define A^i^' := A„ and if A^'''*'^-' 
has been defined, we use aI'''' to denote a i/-dimensional set of the above form 
which contains the maximum number of good sites. Note that if ?/ is as in the 
statement, we have 



|Ai'''| 

because the ratio on the left decreases in i^. 



> V (4-42) 
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Next we set rji, :— — 2"^ + 1 and note that 

2?7i/ - 1 = i]u+i and rj > 77^. (4.43) 

Assuming without loss of gcncrahty that > 1 — 2~'^ — otherwise we can 
take A := in the statement of the lemma — we have 77^ > for all 1 < ly < d. 
We will prove by induction the following claim: For each v — 1, . . . there 
exists a set A^'^^ C GK,ei^) H aI"^' such that 

|A(''^|>r?.|A(r)| (4.44) 

and 

x^yeA'^"^ \xiy,io)-x{x,io)\<2iy[K + e{2n+l)]. (4.45) 

For u = d this clearly implies the desired claim. 

For = 1, we define A^^^ := Ck.A'^) ^ ^" ^- "Hi < this obeys (4.44). 
The bound (4.45) is then a direct consequence of the definition of a good line. 
Suppose now that the claim holds for v and let us prove it for v -\- 1. To this 
extent, let 11 denote the natural projection of A^'^^^' onto A,^i'''' and, given the 
set A''^^ with the above properties, let 

^(''+1) :={a;eGif,e(a;)nAl''+i):n(x)GA('')}. (4.46) 

We now verify that this obeys (4.44-4.45). As to (4.44), the same bound 

for A'''^^ tells us that at most (1 — 77,^)1 A^'^^"'^' | sites in Kn^'^^ do not project 
into A'^"^ . Hence 

> \GKA^)nAl^+'^\-ii~vu)\Ai''-'^\ 

(4.47) 

where we used (4.42) to get the second inequality. In light of (4.43) this implies 
(4.44). 

To prove also (4.45), we pick two sites x,y £ A'-'^^^^ and let x — Il{x) and 
y = n(j/). The claim for v then implies 

\x{y,u;)-xix,uj)\ < 2zy[i^ + e(2n+l)] (4.48) 

while the fact that a; is a good site yields 

I x{x, uj) ~ x{x,uj) \ < K + e{2n + 1) (4.49) 

and similarly for the pair y and y. Combining these bounds and using the triangle 
inequality then implies (4.45) for x and y — with, of course, 21/ replaced by 
2{u+l). □ 

Lemma 4.16 now implies that the corrector is sublinear on average: 
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Proof of Proposition ^.15. Suppose without loss of generality that 5 < S^**, fix 
e < -^5 and note that we can choose K so large that P(0 G G'k,e) > 1 — ^2- 
By the Spatial Ergodic Theorem and ergodicity of P we thus have 

> 1 _ , (4.50) 

once n > uq for some a.s. finite uq = hq^lj). We will assume that uq is so large 
that also 

Sn > md[K + e{2n + 1)] (4.51) 

holds for all n> no- 

By Lemma 4.16, for each n > uq there exists An = An{w) C A„ with 

> (1-2'*,5)|A„| (4.52) 

and (4.40) valid for ah x,y e An. As S < 8~'', we have |A„| > (1-4~'')|A„| while 
\Mn \ A-2,n\ < 4~'*|A2„| = 2"''|A„|. In particular, An n ^ for each n > uq. 
Let ko be the smallest integer such that 2*^° > no and let us pick a site Xk G 
A2k n A2k+i for each k > ko. The bounds (4.40) and (4.51) then give us 

fc-i S '^^^ 

\x{xk,uj)-x{xk„uj)\ < "^\x{xi+i,uj)-x{xe,uj)\ < - ^ 2^+^ < (52'=-2 
e=ko 1=0 

(4.53) 

Choosing fci = ki(uj) > ko so that \x{xko,^^)\ < <>2''^ ^, this and (4.40) imply 

X e A2k ^ \x{x,uj)\ < S2'', k>ki. (4.54) 
But this means that for n € {2'° : k > ki}, 

E l{lx(x,<.)|>M < |An \ < 62''\An\. (4.55) 

|x|<n 

As S was arbitrary, this proves (4.38) for n increasing along powers of two. A 
moment's thought now reveals that the same then holds for the unrestricted 
limit as well. □ 

As for Theorem 4.14, we refer the reader to Biskup and Prescott [18]. It should 
be emphasized that, although the assumptions to all the above are those of the 
annealed invariance principle, we in addition require the validity of the diffusive 
bounds (4.34-4.35). These are by no means guaranteed for a general ergodic P, 
so the problem whether the annealed and quenched invariance principle hold 
simultaneously remains open. 

We close this subsection with a simple exercise concerning the invariance 
principle for the variable-speed continuous-time version of our random walk. 

Exercise 4.17. Suppose the "usual assumptions" and assume that (Xn) obeys 
the Quenched Invariance Principle with the limiting Brownian motion having 
covariance (3.36). Show that the variable-speed continuous time walk Xt obeys 
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a Quenched Invariance Prineiple with the limiting Brownian motion having co- 
variance 

E[{\-Btf] =E(^c^o,4A.«'(c^,a;)]'). (4.56) 

X 

Note that the quantity on the right-hand side is closely related to the infimum 
of (3.9), which was used to define the corrector. The appearance of expectation E 
instead of i?Q is due to the fact that P is invariant for the point of view of the 
particle induced by the VSRW. As to the constant-speed walk, here the quenched 
invariance prineiple follows from the discrete-time case by a strong asymptotic 
concentration of a sum of i.i.d. exponential times. 

4.4- Known results and open problems 

The following sums up the principal steps in the progress towards proving 
quenched invariance prineiple in the class of Random Conductance Models: 

• Strongly elliptic, ergodic P: proved by Sidoravieus and Sznitman [124]. 

• Nearest-neighbor, i.i.d. conductances in d > 2 subject to the conditions: 

E(wo,c) < oo and P(we > 0) > Pc{d) (4.57) 

where Pc{d) is the bond-percolation threshold. Here the quenched CLT has 
been proved in a sequence of papers by Sidoravieus and Sznitman [124], 
Berger and Biskup [11], Mathieu and Piatnitski [101], Mathieu [100], Biskup 
and Prescott [18], Barlow and Deuschel [7] with all approaches synthesized 
together by Andres, Barlow, Deuschel and Hambly [3]. 

• Nearest-neighbor, ergodic P with in d ~ 1,2 with 

E(iuo,c) < oo and E(l/wo,6) < oo (4.58) 

The d = 2 case is proved in these notes; the d = 1 is the content of 
Exercises 3.12 and 4.4 and goes back at least to Kawazu and Kesten [84]. 

We remark that that the condition E(wo,e) < oo is essentially necessary; indeed 
Barlow and Cerny [6] {d > 3) and Cerny [30] {d = 2) proved that for i.i.d. 
nearest neighbor conductances with a-stable upper tail, a < 1, the law of Xnt 
is under proper sealing described by , where Bt is a Brownian motion and 
Wt is the inverse of an independent stable subordinator with index a. In other 
words, the paths are still Brownian but the heavy edges introduce a non-trivial 
trapping effects thus rendering the time parametrization non-linear and, in fact, 
stochastic. We remark that in physics, the limiting process is referred by as the 
fractional kinetics process. 

An important open problem concerns the rate of convergence and quantifi- 
cation of errors in martingale approximations. Although optimal results are 
probably far from reach, interesting ideas have been developed and quantitative 
results derived by Mourrat [106] and Gloria and Mourrat [66]. The aforemen- 
tioned work of Gloria and Otto [67] gives integrability estimates on the corrector 
in d > 3 under strong ellipticity. 
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The Random Conductance Model has been also studied over other base- 
graphs than just Z''. For instance, Caputo, Faggionato and Prescott [26] have 
investigated the random walks over various point processes in M.'^. Independent 
studies for random walks on Voronoi/Delaunay triangulations have been an- 
nounced by Buckley [24] . Ferrari, Grisi and Groisman [58] have constructed the 
harmonic coordinates on such triangulations by means of an interacting particle 
system; namely, a harness process, which is basically a full-space stochastic ver- 
sion of the algorithm described for the finite boxes in Sect. 3.1. The methods of 
Kipnis and Varadhan can be applied even to some deterministic quasiperiodic 
structures; see, e.g.. Teles [128] who recently established an annealed invariancc 
principle for the simple random walk on Penrose tilings. 

Although we are able to control the corrector to the level required for the 
quenched invariance principle, the object itself remains rather mysterious and 
many open questions remain. For instance, regardless of what has been said at 
the end of Section 3.2, the following problem remains of great interest both from 
the perspective of probability and analysis: 

Problem 4.18. Is it true that a.e. realization of random conductances satisfying 
the "usual conditions" admits no non- constant, sublinear harmonic functions? 

Recently, Benjamini, Duminil-Copin, Kozma and Yadin [9] have shown that 
that on the supercritical percolation cluster in Z'^, the space of linear harmonic 
functions is exactly d + 1-dimensional. In particular, a typical supercritical per- 
colation cluster supports no non-constant sublinear harmonic functions. We ex- 
pect this to hold for all i.i.d. nearest-neighbor Random Conductance Models; 
for general environments the problem remains open. 

Another open question concerns the scaling limit of the corrector: 

Problem 4.19. Assume i.i.d. nearest-neighbor, uniformly elliptic Random Con- 
ductance Model. Show that the corrector scales to a Gaussian Free Field. More 
precisely, for any function f : R'^ — >■ K which is smooth and of compact support 
and satisfies J f{x)dx = 0, the law of 



scales, as e I. 0, to a Gaussian with mean zero and variance proportional to 



Progress in the uniformly elliptic case has been achieved in recent work of 
Gloria and Otto [67] who have been able to prove that the corrector is in L'' (P) 
for all g < oo, and thus a tight random variable, in all dimensions d > 3. This 
settled an open problem from [11]. 

Another, perhaps somewhat related, question is that of the very definition of 
the corrector. Indeed, the corrector is defined almost surely for every ergodic law 
on environments P. However, as different ergodic laws are singular with respect 
to one another, it is not clear how to mesh the various correctors together. And 
yet it seems this should be possible: 




(4.59) 



(/,-A-i 
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Problem 4.20 (Universal corrector). Consider the set of nearest-neighbor envi- 
ronments il :~ [a, 6]®^^ where < a < b < oo. Define a function f7 x Z'^ — > 
R'^ such that, for every ergodic law P on Vl, it agrees with the corrector corre- 
sponding to measure P. 

Wc remark that this would be solved if one could find a sequence of local 
functions ifn such that V(/9„ — x almost surely for every P. Note that, although 
may find functions (y9„ for which the convergence takes place in L^ov fo'" 
given P. almost sure convergence requires reduction to subsequences which may 
be strongly P-depcndent. 

The understanding of the Markov chain permits one to consider more compli- 
cated questions. One such question concerns the typical number of points visited 
by the random walk in a given time. This was recently addressed by Rau[118]. 
Another question is the Law of Iterated Logarithm; this was established by 
Duminil-Copin [50]. Next is the question of the behavior of the random walk 
on very thin percolation clusters. This can be studied directly in the case when 
p = Pc where, technically speaking, the percolation cluster does not exist but 
one can still enforce it by conditioning. For the resulting incipient infinite clus- 
ter (lie) in sufficiently high dimensions, Nachmias and Kozma [108] proved 
the Alexander- Orbach conjecture in all dimensions d > 7 — modulo caveats 
regarding the existing level of lace-expansion technology. This conjecture, due 
to Alexander and Orbach [2], states that, on IIC, 

P^"(0, 0) X n^oo. (4.60) 

Notably, this is expected to be false in low spatial dimensions. Related to this 
would be the decay of the diffusive constant for the simple random walk on the 
supercritical cluster for parameter p, as p pc- Here we pose: 

Problem 4.21. Suppose d > 7 and let D{p) denote the limiting variance of the 
simple random walk on the supercritical percolation cluster on TI^ for parame- 
ter p > Pc{d). Show that 

D{p)-{p-p,)\ Pip,. (4.61) 

This problem is closely related to the existence of effective conductivity which 
was studied in, e.g., Grimmett and Kesten [74], Chayes and Chayes [31] and 
Kesten's monograph [85] on percolation. See also Sect. 6. A rather convincing 
argument can be obtained for this by analyzing the formula (3.36) and making 
plausible assumptions on the structural properties of the percolation cluster. 
Resorting to the electrostatic interpretation, the electric current should be car- 
ried only by the backbone of the cluster — which, in the limit p i Pc, becomes a 
"net" of fractal curves. The exponent in (4.61) then comes from realizing that 
in d > 7, these fractals have Hausdorff dimension 2 (although the relation is not 
so straightforward as a simple equality of these numbers). This intuition seems 
be confirmed by observations made in the physics literature, see, e.g., Schr0der 
and Dyre [120]. A main puzzle that remains is whether, and how exactly should 
the exponent 2 in (4.61) be related to the exponent 4/3 in (4.60). 
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We remark that the amount of physics hterature written on this and related 
subject is absolutely overwhelming; just see the articles citing the review by 
Dyre and Schr0der [52]. 

Another very interesting class of applications of the above techniques is the 
random walk in dynamical (albeit still reversible) random environments. We will 
not go into details here, but let us just say that much of Kipnis-Varadhan theory 
carries to this case and so annealed limit theorems are available. However, the 
understanding of quenched invariance principles is far less evolved. Much can be 
said when the dynamics of the environment is Markovian and there is enough 
mixing; one can then get enough control via regeneration arguments. However, 
even here it is far from clear how to formulate convenient, and very general, 
conditions under which invariance principles can be obtained. 

From the perspective of this text, one specific class of dynamical random 
environments is of special interest. Consider a function K — >■ K which is 
twice continuously differentiable and define a collection of coupled diffusions 
i4>x(t))^^z-' via 



where Bt{x) are independent standard Brownian motions. As it turns out, any 
gradient Gibbs measure for the potential V is stationary under this dynamics. 
Assuming that V is convex, and thus V" > 0, we can now define a random 
walk X = (Xt) which at time t at position Xt = x takes a jump to a neighbor y 
at rate V"{(t)y{t) - (j)x{t)). 

An attractive feature of this setting is that it permits us to analyze gradient 
models with convex interactions. For instance, we have the following formula 



for the covariance of the (static) field in two locations with respect to a gradient 
Gibbs measure ^ by means of the expected number of visit to x by the above 
random walk started at — we expect this to be finite only in d > 3 but other 
formulas exists in d = 1,2. Obviously, this generalizes the well-known formula 
from the Gaussian case which is distinguished by the fact that the random walk 
is not coupled to the evolution of the fields. 

The formula (4.63) is one instance of the Helffer-Sjostrand random walk repre- 
sentation of correlation functions for the gradient model. These have been indis- 
pensable in the study of gradient models with convex interactions (e.g., Naddaf 
and Spencer [107], Giacomin, OUa and Spohn [65], Funaki [55], etc). 

5. Heat-kernel decay and failures thereof 

As discussed at length in the previous section, our current strategy of the proof 
of the quenched invariance principle seems to generally require the use of rather 



(4.62) 



y: \y~x\ = l 




(4.63) 
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precise estimates on the probability that the Markov chain moves from x to y 
in n steps. We emphasize that this is conceptually flawed because we seem to 
need a local- CLT type of result to finish a plain CLT. Notwithstanding, the 
study of the heat kernel is interesting in its own right. We will only review the 
techniques that are ultimately relevant for the applications at hand and refer 
to, e.g., the upcoming textbook by Kumagai [90] for a more in-depth treatment 
of that well-developed area. 



5.1. Some general observations 



To set the vocabulary straight, let us first remark that by the heat kernel one 
usually means the quantity 

q„(x,y):=^^^flf . (5.1) 

As one can expect, Pui{x, •) will for large n approach (a multiple of) the station- 
ary measure n^. So q„, being in fact the Radon-Nikodym derivative of P"(a;, ■) 
with respect to Tr^^, is a very natural object to consider. Note that reversibility 
implies q„(x, y) = q„(?;, x). 

Theorem 4.14 required in (4.35) that the return probability generally decays 
as 71^'^/^. It turns out that, should the CLT hold, we cannot hope for a faster 
decay than this: 

Lemma 5.1. Suppose (Xn) is satisfies a CLT with non- degenerate diffusion 
Constanta^. Assume that ir* :— awp^Tr^(x) < oo. Then there is c = c{d,a'^ ,tt*) > 
such that for n sufficiently large, 

Pl"iO,0)>^^UO)- (5.2) 
Proof. We use reversibility and simple estimates to get 

p^"(o,o) = ^pi;(o,x)p:;(x,o) 



Mx) (5.3) 



\x\<^/n 

The sum on the right-hand can be further bounded using the Cauchy-Schwarz 
inequality: 

P?(o.o) > ^^iS(^J|^ («) 

TT* \{x: \x\ < 

But the CLT ensures that P°(|X„| < ^) > ^P{\Bt\ < l/a) for n large, where 
Bt is the standard d-dimcnsional Brownian motion, and |{a;: < v^ll — 

Jqj. gQj-^Q c' ^ c' {d) < CO. □ 
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We remark that a general method of getting such (inchiding "near-diagonal" ) 
lower bounds in elliptic random environments has been put forward by Nash [109] 
and Fabes and Stroock [57]. 

For reasons discussed earlier, the main technical problem is to find natural 
conditions on the Markov chain so that an n^'^l'^ upper bound can be guar- 
anteed. This problem has been studied for over half a century, starting from 
proofs of regularity of elliptic PDEs with irregular coefficients (De Giorgi [40] , 
Nash [109], Aronson [4]) and validity and consequences of Faber-Krahn, Sobolev 
and Nash inequalities for diffusions on manifolds and Markov chains (e.g., 
Varopoulos [130], Carlen, Kusoka and Stroock [28]). A method to get off-diagonal 
bounds — i.e., for q„(a;, y) with x ^ y — has been put forward by Davies [39] 
based on the Came- Varopoulos bound (Carne [29], Varopoulos [129]). 

In the course of time it has been realized that there is a close connection 
between the desired upper bound and the geometric properties of the underly- 
ing state-space. The key property to check is the validity of the isoperimetric 
inequality (Cheeger [32]) or, more generally, the character of the isopertimetric 
profile (Grigoryan [70]). This connection was later transferred to the context 
of (discrete-space) Markov chains by Lawer and Sokal [94] and Jerrum and 
Sinclair [82] (invoking isoperimetric inquality) and, later, by Lovasz and Kan- 
nan [97] and Morris and Peres [99] (based on isoperimetric profile). 

We will not try to delve deeper into the details of historical developments of 
the subject; instead, the reader should consult the many texts that have been 
written on this (e.g., by Coulhon and Grigor'yan [36], Davies [38], Kumagai [90], 
Montenegro and Tetali [105], Varopoulos [131], Varopoulos, Saloff-Goste and 
Coulhon [131], Woess [135], etc). For us the key fact is that with many Markov 
chains we may associate a natural graph structure — simply put an edge between 
any two states in the state spaces that have a positive transition probability of 
a jump from one to the other. This permits us to connect the mixing properties 
of the chain with facts about geometry of this graph. 

To illustrate this on an example, consider a graph that consists of two bulky 
components connected only by a few edges. Clearly, it will take quite a long 
time to exit one component and discover the other. Naturally, one is thus lead 
to comparing the size of a set with the size of its boundary which is expressed 
very well in terms of aforementioned isoperimetric inequalities. 

In what follows we will rely on a result from a recent work by Morris and 
Peres [99] which we find particularly attractive for its probabilistic flavor. Con- 
sider a countable state Markov chain with state space V, transition kernel P 
and a stationary reversible measure tt. For a finite set A G we will measure 
the boundary via 




(5.5) 



and the volume via 




(5.6) 
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Define the function 

'/'M:=mf|^^^^:^(A)<r| (5.7) 

tliat expresses tfie size of the least possible surface-to-volunie ratio for all sets 
with volume less than r. We can call this function the isoperimetric profile. Its 
computation is often facilitated by the following fact: 

Exercise 5.2. Show that in (5.7) we can restrict to A that are connected — in 
the sense that for every x,y Cz A there is a time n with P"(.t, y) > 0. 

We now quote verbatim Theorem 2 of [99]: 

Theorem 5.3. Suppose that P{x,x) > 7 for some 7 G {0,^/2)- For all e > 0, 
all x,y €z V and all n satisfying 



1-7 



4/e 



">1+ -) / -TTTJ (5-8) 

^ 7 ^ J4[,r(x)A7r(.y)] r(t>{ry 

we have 

P^{x,y)<e^{y). (5.9) 

The restriction to uniformly positive holding probability, P(x,x) > 7, is a 
technical nuisance in applications that often requires analyzing a modified chain 
that has this property. 



5.2. Heat kernel on supercritical percolation cluster 

It is quite instructive to check how Theorem 5.3 implies the usual bound for the 
simple random walk and/or elliptic nearest- neighbor environments. However, we 
will instead do something far less trivial; namely, we will show how this theorem 
applies in the case of the random walk on the supercritical percolation cluster. 

Theorem 5.4. Suppose d>2 and p > Pc{d). There is a constant c = c{d,p) < 
00 and a random variable Uq = no(a;) such that for almost every sample of the 
bond-percolation cluster "^^a containing the origin, we have 

For a finite set A C '^oo('^), let d'^A denote the set of open edges in uj with 
exactly one endpoint in A. A simple observation yields 

QiAAn 1 \d-A\ 

7t{A) - 2d \A\ ^^-^^^ 

li u! := 1 for all edges, then d'^A ~ dA. In such circumstances, one has the 
isoperimetric inequality of the form: There is a constant c = c{d) > 0, such that 



\dA\>c\A\^, AdU^ finite. 



(5.12) 
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This inequality cannot hold on because the infinite component contains 
arbitrarily long one dimensional (and other) pieces. However, we can have this 
for connected sets that are not too small compared to their distance to the origin: 

Lemma 5.5. For all d > 2 andp > Pc(d), there are positive and finite constants 
ci = ci{d,p) and C2 = C2{d,p) and an a.s. finite random variable Rq — Roiuj) 
such that for each R > Rq and each uj-connected A satisfying 

A C "^oo n and |A| > (cilogi?)3^ (5.13) 

we have 

\d'^A\>c2\A\^. (5.14) 

There have been a number of proofs of this and/or related results, see e.g. 
Benjamini and Mossel [10], Heicklen and Hoffman [78], Mathieu and Remy [102], 
Barlow [5], Berger-Biskup-Hoffman-Kozma [12], Pete [115]. We will not prove 
this claim here for all p > Pc{d) as the proof uses non-trivial facts from perco- 
lation theory. However, for p very close to 1 there is a much simpler argument 
due to Benjamini and Mossel: 

Exercise 5.6. Show that once p is sufficiently close to one, there is a con- 
stant ci e (0, oo) and a random variable Rq = Ro{ui) < oo such that for all 
R > Ro, 

Ac[-R,R]''nZ'' and |^| > (ci logi?)^^ imply {d'^ A\ > ^\dA\. (5.15) 

Note that from here we will immediately have (5.14) via (5.12). In order to 
see how (5.14) feeds into Theorem 5.3, note that the Markov chain by time 2n 
will not leave the box [— 2n, 2n]. Thus set R := 2n + 1, pick 9 G (0, V2) and for 
A G [—R, R]'^ n Z'* connected let us estimate the ratio in the definition of (p{r) 
by c|^|-i/'' when |yl| > R'^ and by ci?"^ when |^| < R^. (In the second step 
we used that \d'^A\ > 1.) It follows that 

(j){r)>c{r~^/'^ /\R-'') (5.16) 

Plugging this into (5.8), the integral is at most cR"^^ log R + ce~^^'^. This will 
be less than n for e := cn"'^/'^. The inequality (5.10) then follows by applying 
Theorem 5.3. 

A natural consequence of Theorem 5.4 is the result that was first proved by 
Grimmett, Kesten and Zhang [75] by rather different methods (see also Prob- 
lem 1.16): 

Corollary 5.7. The simple random walk on (a.e. realization of) the super- 
critical percolation cluster is recurrent in dimension d = 2 and transient in 
dimensions d > 3. 

Proof. As explained in Sect. 1.3, it suffices to resolve the d = 3 case, but we 
can cover all cZ's just as well. From Lemma 5.1 and Theorem 5.4 we know that 
P[J,(0, 0) X n^'^l'^ . This is summable in dimensions d > 3 and non-summable in 
d = 1,2. The summability is then equivalent to the finiteness of the full-lattice 
Green's function which via (1.41) is then equivalent to transience. □ 
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5.3. Anomalous decay 

From the perspective of nearest- neighbor random walks on Z'', the case of the 
supercritical percolation is a prototype of a non-elliptic situation. However, when 
we think of this walk as the simple random walk on the graph '^^a , it is as elliptic 
as the one can ever hope for. Indeed, any edge in the graph "^^oo has conductance 
one and the ellipticity contrast — the difference between a maximal and minimal 
possible value of the conductance over each edge — is zero. The difficulties in 
the understanding of this walk on "^oo is thus not the lack of ellipticity but the 
intricacies of its random geometry. 

From this point of view it is natural to ask what happens when ellipticity 
gets violated in a robust way. This naturally leads to consideration of i.i.d. 
nearest-neighbor environments where the law of the individual conductances 
is unbounded either from zero or from infinity (or both). The point is that 
both situations can lead to trapping effects although each of them for a slightly 
different reason. We will henceforth focus on the former case and refer to Barlow 
and Deuschel [7] for the latter. 

Suppose, from now on, that the w's are nearest-neighbor, i.i.d. with P(0 < 
Wfc < 1) = 1 and 

essinf(wft) = 0. (5-17) 

Our assumption implies that P(wh G •) has no atom at zero. Thus all nearest- 
neighbor jumps on Z'^ are allowed for the random walk, but some of them may 
be very unlikely. 

It is easy to check that for i.i.d. distribution with these properties, the 
isoperimetry methods sketched above yield a vacuous conclusion. The situa- 
tion becomes even more suspicious after an inspection of the work of Fontes 
and Mathieu [59] in which they design a family of models — not with i.i.d. 
conductances but close enough — in which the expected diagonal heat kernels, 
EP^"(0, 0), decay arbitrarily slowly with n. Of course, this could be just a result 
of taking an average over the environment (remember that we are talking about 
events whose probabilities decay to zero) so one is naturally intrigued by what 
the typical (quenched) decay of Pj"(0,0) might be. 

It will not be too surprising that in d = 1 the trapping can be quite severe 
even for typical uj. Indeed, the following is an interesting exercise: 

Exercise 5.8. Suppose d = 1 and nearest-neighbor, i.i.d. conductances with 
values in (0, 1]. For each sequence A„ — > oo construct a law P such that 

P^"(0,0)>^ (5.18) 

for n large, along a deterministic subsequence Uk ^ oo. 

A moment's thought — and a right idea — then shows that interesting new 
behavior may actually occur even in high-enough dimensions. Consider the fol- 
lowing example from the paper of Berger, Biskup, Hoffman and Kozma [12]: 
Fix a sequence A„ — > oo and define a trap of order n to be the configuration in 
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Fig 5.1. A trap capable of capturing the walk for times ofordern. Here an edge of conductance 
1 is separated by edges of conductance Vn from a path of edges of conductance 1 to the origin. 
Once the walk enters the trap, incurring a cost i/„ of probability, it will stay bouncing back 
and forth for time of order n with a uniformly positive probability. The exit at a particular 
time costs again ^1^. 



Fig. 5.1 and suppose that the distance of this trap to the origin is In. If such a 
trap occurs, we can estimate Pj"(0, 0) as follows. The cost of getting to vertex x 
from the origin is exponential in the distance, i.e., e*^^^"^. Entering the trap at 
the next step costs order V™ of probability. The walk can then be made to stay 
there for the time In minus twice the distance to the trap; this costs only e*^^^^ 
of probability. Exiting the trap at the required time costs one more Vn and the 
trip back to the origin e*^^^"-*. In total, we thus have 

P!,"(0, 0) > eO(^") - eO(i) - eO(^") = '—^ (5.19) 

Now one just beefs up the lower tail of P so that, along a deterministic subse- 
quence rife DO, we have ~ o(logA„). We have a proof of: 

Theorem 5.9. Suppose d > b. For each A„ — )■ oo there exists an i.i.d. conduc- 
tance law P satisfying P(0 < cj;, < 1) = 1, a deterministic sequence nj, — >■ oo and 
a P-a.s. positive random variable C — C{uj) > such that for each n G {nfc}, 

P^(0,0)>^. (5.20) 

n An 

Notice that the above argument yields a similar bound in all dimensions 
d > 2, but this bound has no significant value in dimensions d — 2,3,4 as (by 
the CLT proved by Mathicu [100] and, independently, Biskup and Prcscott [18]) 
Pj"(0,0) decays at least as n~'^/^; cf Lemma 5.1. But in d > 5 this shows that 
the heat kernel may decay more slowly than n~'^^^ and, in particular, there is 
no way that a diffusive heat kernel upper bounds would generally hold. 

An interesting question is whether (5.20) is the worst one can do. The answer 
turns out to be, more or less, in the affirmative: 
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Theorem 5.10. For any nearest-neighbor, i.i.d. conductance law P with P(0 < 
^o,c < 1) = 1 there is a random variable C = C{lu) < oo such that 

d = 2,3, 

P^,"(0,0)<CH <(n-2logn, d = 4, (5.21) 

d>5. 

In addition, we have 

lim n^P^"(0, 0) = 0, P-a.s. in d > 5. (5.22) 

n— >cx) 

and 

lim p2"(0, 0) = 0, P-a.s. in d = 4. (5.23) 

n-J-oo log n 

All except (5.23) in this result is due to Berger, Biskup, Hoffman and Kozma 
[12]; the property (5.23) was derived only recently in Biskup, Louidor, Rozinov 
and Vandenberg-Rodes [16]. The latter group has also shown that, in many cases 
where the heat kernel decays subdifFusively, the trapping phenomenon described 
in the example above actually occurs: the path spends n — o(n) of time in a very 
small spatial region. 

Notice that (5.20) and (5.22) nicely complement each other: anything up to, 
but no worse than, o(n~^) decay can occur in d > 5. A question remains whether 
the log n factor in d = 4 is an artifact of the proof or a real phenomenon. This was 
solved recently by Biskup and Boukhadra [15] who constructed an environment, 
for each sequence A„ — ^ cx), such that 

P^"(0,0)>^, (5.24) 

eventually, along a deterministic subsequence — >■ cxj. The construction is quite 
involved because in d = 4 the trapping occurs more or less equally likely over a 
whole range of exponentially- growing spatial scales (hence the logn factor). 



5.4- Conclusions 



The upshot of the above results and derivations is that with the random con- 
ductance models we are finding ourselves in a somewhat unusual situation when 
the path distribution satisfies a non-degenerate functional CLT and yet the heat 
kernel decays anomalously; i.e., we have a CLT without local CLT. Although this 
may contradict intuition, there is nothing wrong about this: a CLT is a state- 
ment about the bulk of the distribution and a local-CLT is a statement about 
the tails. There is no particular reason why these should match one another. 



6. Applications 

In this section we will try to address some aspects of the applications that 
were introduced in the first section of these notes. Specifically, we will discuss 
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homogenization of discrete parabolic (random) problems, scaling limit of asso- 
ciated Green's functions, convergence of random Gaussian gradient models to 
Gaussian Free Field and, finally, applications to electrostatics. 

6.1. Some homogenization theory 

The phrase "homogenization theory" usually refers to a diverse set of methods 
and ideas that address one of the fundamental problems of material science: the 
computation of macroscopic material constants and characteristics (e.g., heat 
or electric conductivity, resistivity, etc) from the microscopic properties. One of 
the typical mathematical issues resolved by homogenization theory concerns dif- 
ferential equations: Although the microscopic quantities evolve according to an 
differential equation with rapidly varying coefficients, properly rescaled macro- 
scopic versions thereof are governed by equations with smooth coefficients. 

We will not go into the subject and history of homogenization theory in any 
further detail; these can be found in the literature, e.g., the monograph by Jikov, 
Kozlov and Oleinik [83]. Instead, we will attempt to demonstrate the conclusions 
on an example of heat conduction. 

Suppose that some material of a rapidly varying internal microscopic internal 
structure — described at the lattice level of spacing e by a configuration of 
conductances uj — is put in a macroscopic temperature profile at time 0. At the 
lattice level, the evolution of the temperature profile with time is described by 
the Cauchy problem 

d 

—u(t,x):==L^u(t,x), t>0, xeZ^, 

dt (6.1) 

u{0,x)^f{x), xeZ"^, 

where is the operator (1.13) (acting only on the x coordinate) that repre- 
sents the microscopic diffusive properties of the material and / is the initial 
temperature profile. Our first question concerns the existence and uniqueness of 
the solution. We note the classical fact: 

Lemma 6.1. Suppose uj is a sample from an ergodic measure P with KtTi^{0) < 
oo. Let Xt denote the variable- speed continuous-time Markov chain on Z'' with 
generator . Pick f : Z"^ — > M bounded. Then 

u{t,x):^El{f{Xt)) (6.2) 

is the unique solution to (6.1) which is bounded in both t and x. 

Proof. By Exercise 2.8 and the general theory expounded in, e.g., Liggett [95], 
the conditions on w guarantee that a stochastic solution to the backward Ko- 
mogorov equations (1.14-1.15) exits and the semigroup for the VSRW is well 
defined. The fact that (6.2) is a solution is then a consequence of a direct cal- 
culation. Indeed, we have 



u{t,x)= J2 Ki^,^)f{z)- 



(6.3) 
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and the boundcdncss of / and finitcncss of tt^ permit us to exchange the sum 
over z with the time-derivative and L^^. Hence, u satisfies (6.1). 

The remaining issue is thus a proof of uniqueness among bounded solutions. 
Let u{t, x) be such a solution and, for < s < consider the random variable 

M,:^u{t- s,X,) (6.4) 

and let ■= cr{Xr: < r < t). Then {Ms,^s}o<s<t is a martingale. Indeed, 
by the Markov property, on the event {Xg = z} we have 

E^iM,+s - M,|^,) = E'^u{t ~s~S,Xs)- u{t - s, z) 

= E^ [u{t -s-S, Xs) - u{t -s-5, z)] (6.5) 
+ [u{t — s — S, z) — u{t — s, z)] . 

This yields 

lim i [E'::iAL+s - A/,1^,)] = Ufi(i - s, z) - - ^) - 0. (6.6) 

almost surely for every s. Integrating over final intervals and applying the 
Bounded Convergence Theorem proves that {Ms, ^s}o<s<t is a martingale. (At 
s = t we apply continuity from the left.) 

The Optional Stopping Theorem then yields 

EZMo = EZMt (6.7) 

which reads 

u{t, x) = EZuit, 0) = E^^fiXt) = uit, x). (6.8) 

The uniqueness is proved as well. □ 

Exercise 6.2. Construct a configuration of nearest-neighbor conductances on 7L 
for which there is a non-zero solution to (6.1) with u{0, •) :~ 0. 

Our next goal is to describe the asymptotic of the solution for the situation 
when / is a macroscopic profile over a lattice of spacing e. Fix a function / : M'* 
M in I/^''°'^(dx) and let u^'^\t,x) denote the unique bounded solution to (6.1) 
with initial data 

u(0, x) := / dz f{ex + ez), xeZ'^. (6.9) 

"'[0,1]'' 

Under diffusive scaling of space and time, we get the quantity 

u,{t,x) ■.= u^''>{te-^,[xe~^\), i > 0, x £ M'^. (6.10) 

Theorem 6.3. Suppose /: M'' — > M obeys ||/||^2(Rd) + ||V/j|^2(Rd) < oo and 
let V be a law on the conductances satisfying the "usual conditions," (3.23) and 
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(4.3). Let Q denote the generator of the (annealed) limiting Brownian motion 
and let u he the solution to the Cauchy problem 

d 

ttuU, x) = Quit, x), t>0,xe R'^, 

at ' ' ' ' - (6.11) 

Then for each t > 0, 

u,{t,-) — > u(t,-) in L^{dx)(l) L^(F). (6.12) 

Proof. Let Bt be the Brownian motion with generator Q. Then 

uit,x)^E''{f{x + Bt)). (6.13) 
Similarly, resolving the above scaling relations yields 

u,{t, x) := 4^^"J ( / dz f{eXt,-. + ez)) . (6.14) 

By translation-invariance of P and the Cauchy-Schwarz inequality, 

E / \u^{t,x) — u(t,x)\^dx 
Jr 

<f dz f^(\E"^f{x,iz) + eXt,-2)~E'>fix + Bt)\^)dx, (6.15) 

JlO,!]'' J ^ ' 

where x^^iz) := elxe^^\ + ez. 

Our first step is to replace Xf{z) by x in the argument of the first / on the 
right-hand side. The difference tends to zero when e 4- because we have 

J \E^^f{x,iz) + eXt,-2) - Elf{x + eX,,-2)|'dx < 4e2|| V/||i.(R.). (6.16) 

To control the remaining difference, we note that, by the Annealed CLT (in 
analogy with Corollary 4.6) there exists a coupling Q° of the random walk 
and the Brownian motion Bt such that, for any 5 > Q and any t > 0, 

a,{S):^¥.Ql{\eXt,^.-Bt\>5) -^0. (6.17) 
Picking an arbitrary (5 > 0, the bound 

E^QO J\f{x + eXt,-2)^ fix + Btfdx 

<^a,mf\\l,^j,,^+S^Vf\\l..^^,^ (6.18) 
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then shows that the expectation on the left tend to zero as e 4- (foUowed by 
(5^0). The proof is then finished by noting that 

E j (|i;°/(x + eXt,-2)-£;V(a; + Bt)|^)da;<LHS of (6.18), (6.19) 

as imphed by using Cauchy-Schwarz one last time. □ 

Theorem 6.3 exemphfies a statement in homogenization theory. Indeed, a 
solution to the parabolie problem with rapidly varying coefheients does behave, 
at a large scale, as a solution to a parabolic problem with constant coefficients. 
As is seen from Exercise 4.17, the coefficients in the equation, namely, the entries 
in the symmetric, positive semi-definite matrix {qij) in 

Q/(-)-E^.a^ (6-20) 

are given by 

g,, :=E(^wo,4e«-*(w,a:))(e, •^'(w,a;))), (6.21) 

X 

where \I'(w, cc) is the harmonic coordinate discussed at length in Section 3. Notice 
that these are characterized by a variational problem 

d 

=inf E(^wo,4A-a; + V,</j(w))') (6.22) 

where A = (Ai, . . . , A^) S R'^ and where i^j: £7 ^ R runs over all local functions. 
This is the same variational problem that defines the corrector. This is the de- 
sired formula that at least in principle allows us to compute material coefficients 
from its microscopic properties. 



6.2. Green's functions and gradient fields 

The arguments in the previous section can be cast in a more symmetric form 
provided we are willing to invoke some functional analysis. Given an opera- 
tor on ^^(Z'') with coefficients ^{x, y) :— {Sx, ^Sy), we can interpret it as an 
operator on L'^{M.'^) by way of 

{f,^'9)LHw-)-^ f i^{lx\,[y\)f{x)g{y)dxdy. (6.23) 

For any / G L'^{R'^) define 

Mx) e'^^'+'fixe). (6.24) 
In this notation, the statement of Theorem 6.3 implies: 
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Corollary 6.4. For any smooth functions f,g: R — >■ R o/ compact support, 

e-2(<7.,e*'"L-/e)^,(jj,) ^ (5,e*«/)L2(R.), in ^^(p). (6.25) 

Proof. Just note that, in the notation of Theorem 6.3, e~^( gc, e**^ "^"^ fe) = 
{g,u^{t,-)) while {g,c*'^f) — {g,u{t,-)). These tend to each other as u^{t,-) 
u{t,-) in L^{dx)(g) L^{P). □ 

Corollary 6.4 supplies the core idea underlying the proof of our next result: 

Theorem 6.5. Consider any ergodic law P on nearest-neighbor elliptic conduc- 
tances and pick any /, g : — >■ R that are smooth and of compact support. In 
d = 1,2 assume in addition that the integral of f and g over equals zero. 
Then 

(5.,(-U)'V.)^.(K.) ^ (5,(-Q)-V)i.(R.), ini^(P). (6.26) 

Proof. (Sketch) We only sketch the main ideas; details for this setting can be 
found in the work of Biskup and Spohn [19]. All inner products will be those 
in LP'{W^) so we will not make this notationally explicit. 

First let us note that both inner products are well defined. Indeed, — is 
self-adjoint and positive semi-definite with empty kernel (in I'^ilJ^)). Moreover, 
it is invertible on all functions of finite support in Z'^ subject to — in dimensions 
d = 1, 2 — the condition of a vanishing total sum. Uniform ellipticity gives us 
the following inequality between norms: 

( /, (-U)-V) < 11/11^ + c{ /, (-A)- V), (6.27) 

where A is a continuum Laplacian and where the passage from discrete to 
continuum Laplacian is due to [19, Lemma 2.2]. As is not hard to check, replac- 
ing / by on the left and using that ||/e||2 = e^||/||2 while (/, (-A)~Ve) = 
(/, (-A)-V), the bound stih holds ah e < 1. The family in (6.26) is thus 
uniformly bounded. 

By the polarization identity, it suffices to prove the claim for g := /. To this 
end we notice the following representation 



(/„(-U)-Ve) = / (/„e*L-/,)di, 

,oo (6.28) 

= / e~\f,,e'^-'^-f,)dt, 
Jo 

where we scaled t by in the second line. By Corollary 6.4, the integrand on 
the right-hand side tends to (/, e**^/) so, ignoring the important issue whether 
we are able to interchange the limit and the integral, we should have 

/>oo 

(/e,(-U)-V.) -t / (/'C*S)dt. (6.29) 

Jr. 
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The right hand side is again bounded by the fact that Q is uniformly ehiptic, 
and it equals the term (/, (— Q)^^/). 

The key technical point of the proof is thus the control of the tails of the 
integral in (6.28). This is a non-trivial problem where we will have to invoke, 
once again, heat-kernel estimates. This is easier in dimensions d > 3 where it 
suffices to invoke the result of Delmotte [41] : 

Pl{x,y)<^, x,yeZ'',t>0, (6.30) 

with some constant c € (0, oo), uniformly in w — subject to the strong-ellipticity 
condition. For / e L^{dx) this yields 

e-'{f.,e''~''-L)<\\f\\l^. (6.31) 

This is uniformly integrable in all dimensions d > 3. 

In dimension d = 1,2 one needs a corresponding bound on the gradient of 
the heat kernel. Such bounds were proved in the annealed setting by Delmotte 
and Deuschel [42]. Sec Corollary 4.3 in [19] for details. □ 

The above conclusions permit a statement on the random Gaussian field 
introduced in Problem 1.22. Indeed, let {(j)^) be a sample from the Gaussian 
measure with zero mean and covariance (— L^)"^, for a collection of nearest- 
neighbor elliptic conductances k. Recall the notation 4>e{f) from (1.42). Then 
we have: 

Corollary 6.6. Suppose f is smooth with compact support and (in d ~ 1,2) of 
zero total integral. As e i 0, the law of (jj^^f) tends to that of a Gaussian with 
mean zero and limiting variance 

Y^r{Uf)) ^ (/,(-Q)-V),.(«.) (6.32) 

in V -probability, where Q is the generator of the limiting Brownian motion. 

Proof. As 0c is Gaussian, it suffices to prove the convergence of the variances, 
i.e., (6.32). This is (6.26) in disguise. □ 

The key point is that since the limit is non-random, the same will be true 
even if the law of the (/)'s is further averaged over k. This permits the main 
conclusion of the paper of Biskup and Spohn [19] which repharse as follows: 

Theorem 6.7 (Scaling to GFF). Suppose V is as in (1.44) with g compactly 
supported in (0, oo). Let fi be a gradient Gibbs measure for the potential V which 
we assume to be ergodic with respect to the translations of and to have zero 
tilt. Then for every f G Dom((— A)^^/^), the law of (j)^{f) tends to a Gaussian 
with mean zero and covariance 

<^}--^{f.{-Q)-'f)m^.y (6.33) 
where Q^^ is the inverse of the operator (6.20). 
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Proof. (Sketch) Conditional on the k's, the law of the i^'s is Gaussian. Hence 
i5^(e'^.(/)) =i5^(i?^(e'^=(/)|.c)) 



iB,,(0,(/)l«)-5Var(0.(/)|K)\ 



(6.34) 



The above tells us that Var(0e(/)|K) — > cr^ in probability, but we cannot use it 
unless we can simultaneously deal with the conditional mean Ef^{(j)^{f)\K). The 
most substantive part of the result is the following representation: Suppose /i 
has tilt t. Then 

E^{q^^ -(j)y\K)=t- [*(k, x) - ^{k, y)] , (6.35) 

where 4" is the harmonic coordinate. This proves that the conditional mean is 
identically zero when t = 0. Then (6.34) gives 

R (ci^=(/)) e-^'/. (6.36) 

Since this holds for all multiples of / as well, Levy's characterization of conver- 
gence in law implies the desired claim. □ 

We note that the above Gaussian field with random (ergodic) covariance 
structure has been (probably first introduced and) studied by Caputo and 
loffe [27, Section 4.5]. Their motivation was to provide a link between the deriva- 
tive of the exponential rate function for changing the tilt of the field — the so 
called surface tension — and the diffusivity of the corresponding random walk 
among random conductances. For the above Gaussian case, this link is veri- 
fied by a direct calculation, but for general uniformly-convex interactions — for 
which one still has a random- walk representation (Naddaf and Spencer [107], 
Giacomin, OUa and Spohn [65]) — it remains conjectural despite serious effort. 



6.3. Random electric networks 



Theorem 6.5 can be understood as an application of homogenization theory to 
electrostatic equilibrium. Indeed, given an assignment of charge p{x) at vertex x, 
we wish to find an electrostatic potential cp^^ : Z'' — > M satisfying the Poisson 
equation 

lu^vM) ^ Pix), xeZ^ (6.37) 
with the normalization V'aj(O) = 0. As to the existence of solutions, we have: 
Lemma 6.8. Suppose p G r)om{{—Li^)^^^^) which is equivalent to 

sup (p, (e - U)~ V)<>2(zd) < oo- (6-38) 

Then {tp: £{(p) < oo} contains exactly one function ipi^ satisfying (6.37) and 
V'w(O) = 0. Moreover, we have 

inf {'?(¥') + {'P,p}- Si-p) < oo} = -isup(p, (e- L^)-"" p) (6.39) 
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and every minimizing sequence ip^"^"^ on the left has the property that £{ip^"''^ — 
^u>) 0. In particular, the infimum is achieved by Lp^^. 

Proof. (Sketch) The containment p e 'Dou\{{—L^)~^''^) guarantees that the 
suprcmum in (6.39) equals (p, {—\-^)^'^p)i2(^idy A completion of square yields 

£{^) + p)^£{p^ {-^.r"^p) - \ (p, (-U)- V),.(^.) (6.40) 

which by the fact that the Dirichlet energy is non-negative implies that > holds 
in (6.39). Let thus (/s^"^ be a minimizing sequence. The parallelogram law (see 
(3.28)) then immediately gives that£:((^(")-(p('")) ^ as m, n — >■ oo. Moreover, 
£'((p'")) must remain bounded because if we had £(iy9'"-') — > oo, then (<^("\/9) 
would tend to — oo at the same rate as £'(y5^"-') tends to +00. But this is not 
possible as, by p € Dom((— L^)~ ''2) and Cauchy-Schwarz, \{ip^'^'^ , p)\ grows at 
most as the square root of £'((/'^"^). 

Passing to m — > 00 we thus construct a minimizcr ^p^ on Z'' with £{'Pui) < 00. 
By adding small perturbations, wc find that (poj solves (6.37). The identity (6.40) 
and the fact that S{p) = only for constants then shows that puj with a 
prescribed value at one lattice site is unique. □ 

Having dismissed the questions of existence and uniqueness, let us now in- 
vestigate what happens when we scale the lattice to have spacing e and scale 
the charge density to maintain a fixed macroscopic profile. As we will see, the 
following is just a rewrite of results proved earlier: 

Theorem 6.9. Suppose uj is a sample from an ergodic measure on elliptic 
nearest-neighbor conductances and let Q denote generator of the (annealed) lim- 
iting Brownian motion for this environment. Suppose f : ^ R is smooth and 
of compact support with J f{x)dx ~ and let 

p,{x):=[ f{ex + €z)dz. (6.41) 

Define 

Vpi;'(^):=e'(LjV)(L^AJ), ^eR". (6.42) 

Then pl^^ — >■ (p weakly in L^(dx) in P -probability, where p is the solution to the 
Poisson equation 

Qf,(x) = f{x), x e R'^. (6.43) 
(This is well defined as f ^ Dom((— Q)""'^). 
Proof Let g: R'' M. Then 

(ff,V'i;^)i.(K.)= E / dze''g{ex + ez)^^^\ex) (6.44) 

where we used that p^'^'>{ex -\- ez) = p'-'^^{ex) for all x G Z'' and all z G [0, l)'^. 
Invoking the definition of and some elementary rewrites, we then get 

(5,<^L'^)i2(R.) = (5„(U)-Ve)i2(R.)- (6.45) 
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By Theorem 6.5, the right-hand side tends to (g, Q ^ f)L^(Rd-j 
thus in probabihty. It follows that ip[^'' — > Q^^f weakly in 



) in 
idx). 



(P) — and 
□ 
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